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DECOMPOSITION OF THREE-DIMENSIONAL STEADY NON-ISENTROPIC 
COMPRESSIBLE EULER SYSTEM AND STABILITY OF SPHERICALLY 
SYMMETRIC SUBSONIC FLOWS AND TRANSONIC SHOCKS UNDER 
MULTIDIMENSIONAL PERTURBATIONS 

LI LIU, GANG XU, AND HAIRONG YUAN 


Abstract. We develop a method that works in general product Riemannian manifold to decom¬ 
pose the three-dimensional steady full compressible Euler system, which is of elliptic-hyperbolic 
composite-mixed type for subsonic flows. The method is applied to show stability of spherically 
symmetric subsonic flows and transonic shocks in space R 3 under multidimensional perturbations 
of boundary conditions. 


Contents 


1. Introduction 

2. A decomposition of steady Euler system 

2.1. Some well-known reductions 

2.2. Second-order equation for pressure 

2.3. The decomposed system and its equivalence to Euler system 

3. Existence and uniqueness of spherically symmetric subsonic flows 

3.1. Formulation of subsonic flow problem for steady compressible Euler system 

3.2. Existence of spherically symmetric subsonic flows 

3.3. Global uniqueness of spherically symmetric subsonic flows 

4. Stability of spherically symmetric subsonic flows 

4.1. Problem (S2) 

4.2. Problem (S3) 

4.3. Proof of Theorem 14.11 

5. Stability of spherically symmetric transonic shocks 

5.1. Transonic shock problem (T) and main result 

5.2. Reformulation of R-H conditions 

5.3. Problem (Tl) 

5.4. Problem (T2) 

5.5. Problem (T3) 

5.6. Problem (T4) 

5.7. A linear second order nonlocal elliptic equation with Venttsel boundary condition 

5.8. Solvability of Problem (T4) 

6. Appendix 


2 

4 

5 

5 

6 

7 

8 

Ji 

10 

12 

12 

15 

17 

19 

20 
21 
23 
26 
27 

29 

30 
30 
42 


Date : May 13, 2015. 

2000 Mathematics Subject Classification. 35M30, 35R35, 35Q31, 76H05, 76L05, 76N10. 

Key words and phrases. Euler system, three-dimensional, elliptic-hyperbolic composite-mixed type, decomposition, 
subsonic flow, transonic shock, stability, uniqueness, Venttsel condition, nonlocal elliptic operator. 


1 







2 


LI LIU, GANG XU, AND HAIRONG YUAN 


6.1. Solvability and estimate of transport equations 

6 . 2 . Solvability and estimate of div-curl system on sphere S 2 

6.3. Estimates of higher-order terms 
References 


42 

43 
43 
46 


1. Introduction 


We are interested in the three-dimensional steady non-isentropic compressible Euler system that 
governing the stationary motion of subsonic flows of perfect gases. Let p, p, and s represent the 
pressure, density of mass, and entropy of the flow respectively. To be specific, we consider polytropic 
gases, namely the constitutive relation is given by p = A{s)p^ , with 7 > 1 being the adiabatic 
exponent, and H(s) = ko exp(s/c u ), for two positive constants ko and c v . The sonic speed is 
c = V'yp/p- Let u Le the velocity of fluid flow, which is a vector field whose integral curves are 
called fluid trajectories. The flow is called subsonic if |u| < c, and supersonic if |it| > c. Now, for 
‘div’ and ‘grad’ being respectively the divergence and gradient operator, set 


di v(pu < 8 > u) + gradp — pb, 

( 1 . 1 ) 

di v{pu), 

( 1 . 2 ) 

di v(pEu) — pb ■ u — pr. 

(1.3) 


Here E = ^\u\ 2 + is the so called Bernoulli constant m p. 22 j), which is the total energy of 
the fluid per unit mass, and b is the exterior body force vector acting on the fluid (per unit mass), 
and r is the heat supply (per unit mass and unit time). Then the steady, three-dimensional full 
Euler system for compressible fluids with exterior force b and heat supply r reads (c/. [13l (3.3.29) 
in p.62]): 


V = 0 , ¥>1 = 0 , cp 2 = 0 . 


(1.4) 


They represent respectively the conservation law of linear momentum, mass and energy. 

The Euler system (11.41) fits into the general form of multidimensional balance laws. For supersonic 
flows, it could be written as a symmetric hyperbolic system m Section 2.1], for which a local well- 
posedness theory of (piecewise) classical solutions is now available (see [2] [13]). For subsonic flows 
without stagnation points (i.e. 0 < |tt| < c), system (11.41) is of elliptic-hyperbolic composite-mixed 
type El p.538]. In fact, using the standard Descartes coordinates (z 1 , z 2 , z 3 ) of Euclidean space 
M 3 , one could write (11.41) (with 6 = 0 and r = 0) in the matrix form 


A 1 {U)d z iU + A 2 {U)d z 2U + A 3 {U)d z zU = 0. 

Here U = (u,p, s) T € M 5 , and u = ( ui,u 2 ,u 3 ) t € M 3 ; for n = (ni,n 2 ,ns) T € M 3 , the 5x5 real 
symmetric matrices A \, A 2 and A 3 are given by 

3 / ^ n T 0 \ 

I pc 

2 _^Aj(U)nj = n p(u ■ n)I 3 0 , 

i =1 y 0 0 u ■ n j 

where I 3 is the 3x3 identity matrix. If u 1 € (0, c) and |w| < c, then for any (£, 77 ) € M 2 \ {(0, 0)}, 
A\{U) is nonsingular and Ai{U)~ 1 {A 2 {U)E, + A 3 (U)rj) is diagonalizable over C, with a pair of 
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conjugate imaginary eigenvalues A± = Xr ± \J—1 A/ of multiplicity one, and a real eigenvalue 
A 0 = U2 ^^ 3TI of multiplicity three. Here 


Ar 



(u 2 ^ + uzrj) and A / 


( c 2 ~ C»i ) 2 )(£ 2 + 7 ? 2 ) ~ (^ 2 ! + u 3 r ]) 2 

(ui ) 2 — c 2 


are real. A system of first-order partial differential equations with such properties is called of 
elliptic-hyperbolic composite-mixed type. To our best knowledge, up to now there is no any theory 
for multi-dimensional systems of such nonclassical type. (See [T2] for some results about the two- 
independent-variable case, which is based on the theory of generalized analytic functions.) 

The basic idea to treat such elliptic-hyperbolic composite-mixed type system is to decompose it 
to an elliptic-hyperbolic coupled system, which consists of several hyperbolic equations coming from 
the real eigenvalues, and some elliptic equations originating from the imaginary eigenvalues. This 
was proposed by S. Chen in [SJ [9]J in studying a transonic shock problem for the two-dimensional 
steady compressible Euler system. He achieved this by introducing Lagrangian coordinates via 
conservation of mass and using characteristic decomposition method. This was then generalized by 
other authors, and paved the way for a quite well developed theory of steady subsonic flows and 
transonic shocks in two dimensional case (see, for example, ed u EH ms nn and references 
therein). 

An effective decomposition for the three-dimensional Euler system (11.41) is apparently more dif¬ 
ficult. The third author proposed a method in [30] to decompose the Euler system with cylindrical 
symmetry. This was further developed in ed, where the authors used characteristic decompo¬ 
sition and pseudo-differential calculus to decompose (11.41) to a second order elliptic equation of 
pressure coupled with four transport equations of velocity and entropy, and solved a problem on 
the instability of a class of transonic shocks in three dimensional straight duct under perturba¬ 
tions of the back pressure (he., the pressure given on the exit of the duct). (The case of three 
dimensional straight duct with general smooth cross section was considered by S. Chen in [10].) 
The merit of this problem is that the background transonic shock solution is piecewise constant, 
which simplifies the decomposition since many terms can be considered directly as higher order 
terms in linearization. However, if we wish to deal with non-constant background solutions, such 
as spherically symmetric subsonic flows or transonic shocks established in ED, which turns out to 
be very important in understanding transonic shock phenomena in divergent nozzles EDGE] , we 
need to refine the decomposition to handle these extra terms. In [7], the authors proposed another 
method to decompose the Euler system (11.41) . However, unlike ED, it is quite difficult to show 
that the decomposed problems are equivalent to the Euler system, since there involve too much 
differentiations of the Euler system. 

In this work, we combine the ideas presented in ED ID to decompose the Euler system (11.41) to 
a second order equation of pressure plus four transport equations of entropy, Bernoulli constant, 
and tangential velocity components respectively, and show that under appropriate boundary con¬ 
ditions, for quite regular solutions, the decomposed system is equivalent to the Euler system (jl.4l) 
(see Theorem ED. Then we apply this method to study two typical problems in mathematical 
gas dynamics, namely the stability and uniqueness of spherically symmetric subsonic flows and 
transonic shocks under multidimensional perturbations of suitably given boundary conditions. The 
main results are Theorem ED and Theorem l5.ll It is remarkable that as in [7j, we utilized some tools 
from differential geometry to avoid the topological singularity of the spherical coordinates of M 3 , 
and our decomposition works in general product Riemannian manifold. In the studies of transonic 
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shocks, we discovered many beautiful intrinsic structures of the Rankine-Hugoniot conditions and 
the Euler system, and there arise problems like solving div-curl system on sphere (c/. (15.1161) 1. and 
a nonclassical nonlocal elliptic equation with Venttsel boundary condition (c/. (I5.77p l. It should 
be noted that the methods and ideas presented in this work can be used to treat some other typical 
gas dynamics problems involving subsonic flows. We also remark that due to the carefully designed 
decompositions, we avoid loss of derivatives and the nonlinear (free boundary) problems could be 
solved just by using a standard generalized Banach fixed point theorem. 

We note that for subsonic flows in a straight duct with rectangular cross section, Chen and 
Xie gj found an approach different from ours to deal with isentropic Euler system under certain 
boundary conditions on the entry and exit. We also recommend the paper [28] of Weng for many 
interesting observations on the equations (11.41) . The papers [11, [10] 0 T] [28j are the only works we 
know investigating the genuinely three dimensional steady subsonic Euler system. 

We remark in passing that, as a byproduct, we also proved in Theorem 13.II the global uniqueness 
of spherically symmetric subsonic flows in the class of C 2 functions; namely, any C 2 solution to the 
Euler system (jl.jp with suitable spherically symmetric boundary data must be itself spherically 
symmetric. It is reduced to the case of irrotational flow by studying the transport equations of 
vorticity. However, similar global uniqueness result for the transonic shocks is unknown, even for 
the simpler piecewise constant case studied in DU. or the two-dimensional cylindrical symmetric 
case (the case studied in [21]). (Note that the global uniqueness of several piecewise constant 
transonic shocks had been proved for two-dimensional steady compressible Euler system in m-) 

The rest of the paper is organized as follows. In Section [2] we prove a general theorem (Theorem 
m on the decomposition of Euler system (11.41) . In Section [3] we formulate the problem of stability 
and uniqueness of spherical subsonic flows, and prove Theorem 13.11 on the global uniqueness. In 
Section g] we apply Theorem 12.11 to prove Theorem 14.11 Section [5] is devoted to formulating the 
transonic shock problems and proving Theorem 15.II on the stability of certain spherically symmet¬ 
ric transonic shock solutions under general perturbations of upcoming supersonic flows and back 
pressure. Finally, Section [6] is an appendix, which contains solvability and estimates of transport 
equations, and a div-curl system on sphere. We also present some details of estimates of higher 
order terms in suitable function spaces. These higher order terms appear during linearization, and 
their estimates turn out to be quite straightforward. 

2. A DECOMPOSITION OF STEADY EULER SYSTEM 

We develop here a strategy to decompose (II.4|) . the three-dimensional steady full compressible 
Euler system, into a second order equation for pressure and four transport equations. Then we 
show that any quite regular solution to the decomposed system is also a solution to the Euler 
system. 

For convenience of applications in mind, we would like to write the decomposition in a quite 
general way by using some terminologies from differential geometry. For a given Riemannian 
manifold A4 with metric tensor G, we use d to denote the exterior differential operator and D u uj 
the covariant derivative of a tensor field cj with respect to a vector field u in AT C u oj is the Lie 
derivative of u> with respect to u. As a convention, repeated Roman indices will be summed up 
for 0,1,2, while repeated Greek indices are to be summed over for 1,2. We also use standard 
notations such as C k , C k,a (Q) and H S (Q) to denote respectively the class of £;-times (k £ NU {0}) 
continuously differentiable functions, and the Holder space of C k functions on an open (or closed) 
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set fl whose all k -th order derivatives are Holder continuous in with exponent a € (0,1], and the 
Sobolev space W s,2 {yt) with s > 0. 


2.1. Some well-known reductions. We firstly derive some well-known equations from the Euler 
system m which are valid only for C 1 flows without vacuum (p > 0). 

The conservation of mass p\ = 0 can be written (equivalently for C 1 flows) as 

ipi = D u p + p div u = 0. 


Then the conservation of energy <p 2 = 0 becomes the Bernoulli law 

+3 - ~(<P 2 ~ E(p i) = D U E - G(b, u)-r = 0. (2.1) 

P 

By the identity div(u <8> v) = (divu)rt + D v u for two vector fields u and v, the conservation of 
momentum ip = 0 turns out to be 

ipo = —Up — ipiu) = D u u -|—gradp — b = 0. (2.2) 

P P 

Since \u\ 2 = G(u,u), and ^D u (ju\ 2 ') = G(D u u,u) = G(<po,u) — j } D u p + G(u,b), we have 

Y (^3 - G((po,u)) = D u A(s ) - ^~T r = 0; (2.3) 

pi 1 ' 7 pi 1 

namely, for r = 0, the well-known fact that the entropy is invariant along flow trajectories for C 1 
flows. From this, the conservation of mass may be written as 


_ a ¥>i , ^4 

Vi — -1-tun 

P 7 Ms) 


divn + 


DuP 

7 P 



(2.4) 


We now derive a transport equation for vorticity. For a vector field u = u l di in A4, we use 
u = u^Gijdx 1 to denote its corresponding 1-form. Then by the formula C u u = D u u + d(-^j-), (12.21) 
is equivalent to 

po = D u u + — - b = C u u - d f + — -6 = 0. (2.5) 

P \ 2 J P 

Noting that d commutes with Lie derivative and d 2 = 0, this implies 


C u (du) = -d 



A dp + d6. 


( 2 . 6 ) 


2.2. Second-order equation for pressure. Straightforward computation yields the following 
tensor identity: 

di v(D u u) — D u (divu) = C\C\{Du <g> Du) + Ric(u, u), (2.7) 

where Du is the covariant differential of the vector held u, Ric(-,-) is the Ricci curvature tensor, 
and Cj(T ) is the contraction on the upper i and lower j indices of a tensor T. From (12.2|) and (|2.4|) . 
we get 

di v(D u u) - D u (divu) = divp 0 - D u (p 1 + D u (^~ - 7 1 r) - div(^—^ - 6). 

7P c z P 
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It follows the identity: 


D u ipi - div^o 


D, 




V IV J \ P J 

—C\C\{Du (8) Du) — Ric(u, u) — (7 — 1 )D U ^ + div6. (2.8) 

We remark that in many applications, M is a connected open subset of M 3 with metric G induced 
from the standard metric of R 3 , so Ric(u, u) = 0. 


2.3. The decomposed system and its equivalence to Euler system. We need some global 
information of M to show that any regular solution of suitably combined reduced equations is also 
a solution of the Euler system (11.411 . To be specific, we assume that M is a product Riemannian 
manifold given by At = (a, b ) x M, where M is a closed surface, and (a, b) C M is an open interval. 
Then 8A4, the boundary of Af, is given by Mq U Mi, with Mq = {a} x M and M\ = {6} x M. We 
also suppose that (a; 0 ,# 1 , x 2 ) with x° € (a, b), {x l ,x 2 ) € M is a local coordinates system on M. 
The following is the first main theorem of this paper. 

Theorem 2.1. In addition to the above assumptions onM, for polytropic gases, suppose also that 
p € C 2 (Ai) n C 1 (M), p,u € and p / 0, u° = G(u, do) / 0 in M. Then p, p, u solve the 

Euler system m in A4 if and only if they satisfy the following equations in A4: 

D u A(s) - ^-r = 0, (2.9) 

D U E — G(b,u) — r = 0, (2.10) 

D u ~ div _ C\Cl{Du <gi Du) - Ric(u,u) - (7 - 1)D U Q + div6 

+L°( : —^ + divu — —r) + L l (D u E — G(b, u) — r) 

7 p c z 

+L 2 (D u A(s) - 2^1 r ) + L\D u u + - b) = 0, (2.11) 

P 7_i P 

G(D u u + ^^--b,d a ) = 0, a = 1,2 (2.12) 

P 

and the boundary condition 

— + divu - 'K^r + L X [D U E - G{b , u) - r) + L 2 {D u A{s) - ±r) 

7 p c p 7 1 

+ LAD u u + — — b) = 0 on Mi. (2.13) 

P 

Here L°(-) is a linear function, and L k {-), Tfc(-) are smooth functions so that L k ( 0) = 0, Tfc(0) = 0 
for k = 1,2,3. 

Proof. The necessity is obvious from the above deductions. For sufficiency, by comparing (12.91) with 
(12.31) , (12.101) with (12.11) , (12.111) with (12.81) , (12.121) with (12.21) , and (12.131) with (12.41) , we conclude that 
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the above equations (12.91) (|2.13D are equivalent to the following expressions: 

(£>4 = 0 = ip 3 - G(tp 0 ,u) in M , (2.14) 

p> 3 = 0 = ip 2 — Eipi in M, (2-15) 

D u <Pi - div¥> 0 + £°(<Al) + L 1 ^) + L 2 (ip 4 ) + L 3 (tp 0 ) = 0 in M, (2.16) 

G(ip 0 ,d a ) = 0, a = 1,2 in M, (2.17) 

(pi + Li(ip 3 ) + L 2 {ipa) + ^ 3 (^ 0 ) = 0 on Mi. (2.18) 


From (12.141) and (12.151) we see ¥>4 = 0, <^ 3 = 0, hence G(ipo,u) = 0. Set u = u k dk■ Then (12.171) 
yields G((po,u°do ) = 0. Note that u° never vanishes in M. by our assumption, one concludes that 
G(<po,do) = 0, hence cpo = 0. Therefore (12.161) is reduced to D u <pi + L°(<pi) = 0, and (12.18[) becomes 
to be if 1 = 0 on M\. Recalling that L° is a linear function, these imply that yd = 0 in A4. So by 
ipi = 0, we have ipi = 0, and hence from <^3 = 0 to get <p 2 = 0. Finally, by (12.21) . it follows that 

= 0. So equations (jl .4|) hold as desired. □ 

Remark 2.1. It is clear that we could propose the condition (|2.13l) on Mo and the same conclusion 
holds. Later in the studies of transonic shocks we will actually restrict (12.131) on a shock-front. All 
the functions L k ,L^ in the theorem are to be specified in the applications. 

Remark 2.2. We note that unlike [7], the regularity of pressure here is one order smoother than 
other unknowns. An iteration scheme works to solve the nonlinear problem, due to the special 
structure of the equation ( 12 . 111 ) . where second order derivatives of pressure occur, while for other 
unknowns (such as p,u), only first order derivatives appear. 

Remark 2.3. We may consider the case that M is a closed manifold, such as T 2 , the flat 2-torus 
given by M 2 /Z 2 (Z is the set of integers), or S 2 , the unit 2-sphere in R 3 . For the former case, the 
existence and stability of constant subsonic flow had been studied in [28j and [3] (for isentropic 
gas), and the transonic shock problem was solved in |llj . These could be handled in a similar 
but simpler way than the latter case, since the special solution is (piecewise) constant. So in this 
paper we focus on the case that M = S 2 , for which the background solution is variable, and many 
interesting new phenomena, like nonlocal elliptic problems, might occur (see Section [5]). 

Remark 2.4. We consider in this section decomposition of (11.41) with exterior force and heat supply 
for convenience to applications to other problems later. In the rest of this paper, like all the 
references mentioned above for subsonic flows and transonic shocks, we only consider the Euler 
system (11,4|) with exterior force 6 = 0 and heat supply r = 0. 

3. Existence and uniqueness of spherically symmetric subsonic flows 

In the following two sections, we apply results in Section 12.11 and Theorem 12.11 to a physically 
interesting problem, namely, the uniqueness and stability of spherically symmetric subsonic flows 
in R 3 . We firstly formulate rigorously the boundary value problem to be concerned, then show 
the existence of spherically symmetric subsonic solutions under suitable boundary conditions, and 
prove that such solutions are unique in the class of C 2 flows. Then in Section [4] we prove that 
some of these spherically symmetric subsonic solutions are stable under general multidimensional 
perturbations of boundary conditions. 
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3.1. Formulation of subsonic flow problem for steady compressible Euler system. We 

now specify A4 = (ro, n) x S 2 for two constants 0 < ro < r\ < oo. So A4 is a spherical shell in M 3 . 
Let r = x° € (ro,ri), and x' = ( x l ,x 2 ) be a (local) spherical coordinates on S 2 . The Euclidean 
metric of M 3 in the coordinates (x°,x') takes the form G = G^dx 1 <g) dx 7 = dx° <g> dx° + (x°) 2 g, 
where g = dx 1 (g> dx 1 + (sinx 1 ) 2 dx 2 <S) dx 2 is the standard metric of S 2 . For later reference, we list 
below all the nonzero Christoffel symbols associated with G (note that F) fc = TJT): 

r° — _ T ° r° — _T^tsin x 1 ! 2 r 1 — r 1 —— 

i n — x , r 22 — x l* 5111 x ) i 1 01 — 1 10 — c o ’ 

r^2 = sinx 1 cosx 1 , ro 2 = r| 0 = -^, r 2 2 = rli = cotx 1 . 

We also set \[G = (x°) 2 sinx 1 , and (G*- 7 ) = (Gjj) -1 . 

Now write the velocity as u = u°do + u! , with v! a x°-dependent vector field on S 2 , which is 

u' = v}d\ + u 2 d 2 in a local coordinates system. We call u° the normal velocity and v! the tangential 

velocity. We prescribe the pressure on Mq: 

p = p 0 (x'), (3.1) 

and the following boundary conditions on M \: 

E = Ei(x'), s = si(x'), u = Ui(x'). (3-2) 

Here u'^x') is a given vector field on S 2 , while E\, si and po are given functions on S 2 . 

Problem (S): Solve the Euler system (11.41) in A4, subjected to the boundary conditions 

m- D- 


Remark 3.1. Choosing suitable boundary conditions for the steady subsonic compressible Euler 
system to formulate a well-posed problem is a delicate issue both for theoretical studies and nu¬ 
merical computations, due to its nature as a nonlinear elliptic-hyperbolic composite-mixed system. 
It had been shown that some choices of boundary conditions, such as given pressure both on Mo 
and Mi, will lead to ill-posed problems, even in the two-space-dimension case (cf. [32]). The 
conditions (13.ID (13.2j) are compatible with the decomposition stated in Theorem 12.II They are also 
motivated by previous studies of transonic shocks and subsonic flows (cf. [32, EU ED] etc.). For 
example, giving a boundary condition of pressure is a physically interesting case in studying nozzle 
flows m p.373, p.385]; while the Bernoulli constant E and the entropy s are usually determined 
by the upstream or downstream far-held flows. 

3.2. Existence of spherically symmetric subsonic flows. We now consider the following spe¬ 
cial case of Problem (S). 

Problem (SI): Solve the Euler system ([1.4)1 in A4, subjected to the boundary conditions 
(EID-O, where po, Ei, si are all constants, and vf = 0. 

We can construct a special symmetric solution to this problem. Suppose that the solution 
depends only on r, and the tangential velocity v! is identically zero in A4. Then the Euler system 
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(11.41) is reduced to the following ordinary differential equations [ 311 : 

d u 2 c 2 u 

dr r(u 2 — c 2 ) ’ 

dp _ 2 pu 2 

dr r(u 2 — c 2 ) ’ 

dp 2 pc 2 u 2 

dr r[u 2 — c 2 ) 


(3.3) 

(3.4) 

(3.5) 


Here, for simplicity, instead of u°(r), we have written u = u(r) to be the normal velocity. Let 
M = u/c be the Mach number. Then it solves 


d M _ lM[2 + ( 7 - 1 )M 2 ] 
dr r M 2 — 1 


So the flow is always subsonic if it is subsonic at the entry r = ro and u is positive in A4. Integrating 
(13.61) yields that 


r 


[ 2 + ( 7 -l )M 2 Y 
Cl y/M 


1 

2 ( 7 - 1 ) 


(3.7) 


Also note that 


d u 2 u 

dM “ 2M + (7 — 1 )M 3 ’ 


hence 


M 

U = Co — _: 

V 2 + (7 -1)^ 2 


Both C 1+2 are constants to be determined by boundary conditions. 


Lemma 3.1. Suppose that 

Po>(--- ! -—j] L , 1 7 _ 1 A(si)T rT , u°(r 0 ) > 0. (3.8) 

\7 7 + 1/ 

T/ien t/iere is one and only one symmetric subsonic solution U to Problem (SI). Furthermore, the 
solution is real analytic. 


Proof. 1. From (12.91) and (12.101) . we see E = E\ and s = si in M. So the assumption (13.81) means 
that the flow is subsonic at the entry r = ro. By po and 7 we could solve p = po > 0 at r = ro 
from p = A(s)p' y , and then u = uq = u°(ro) > 0, M = Mq € (0,1) at r = ro by the fact that 


Ei = lu 2 + 


ip 


2“ 1 (7—1)+ 

2. Now from (|3.7I) . we could determine ci so that 


r = r o [2 + ( 7 -l)M 2 p +2 +-i) . 

[2 + (7 — 1)Mq]^ +5 L^T) y/M 

For r > ro, we have 0 < M < Mq. 

3. Then, we solve (13.31) (13.51) for r € [ro, ri] with initial data p = po,p = po, u = uq at r = ro, to 
get the unique subsonic solution U_(r ) by the theory of ordinary differential equations. Obviously 
the solution is smooth (actually real analytic), and u > 0 in M. □ 
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3.3. Global uniqueness of spherically symmetric subsonic flows. We now prove that the 
spherically symmetric subsonic solution to Problem (SI) constructed above is unique in a larger 
class of C 2 functions. To this end, we need the following lemma on vorticity. 

Lemma 3.2. Suppose that p, p and u = u l di are C 1 and solve the Euler system (ED. with E and 
s being constants in A4. The vorticity of the flow is given by 

du — -LOkmdx ' A dx , LOkm — dfcfGjmU ) 9 m (Gj/ c 'U ), 

where u = u l Gijdx 3 is the 1-form associated with the vector field u. Then u km = —ai m k, and for 
each k, there holds in A4 the algebraic identities: 

u m u km = 0. (3.9) 

Proof. 1. For s a constant in AA, the constitutive relation may be written as p = Ap 1 , with 
A = ko exp(s/cjy) being a constant. 

2. From E = \\u\ 2 + and using c 2 = 7 p/p, we have \d{\u\ 2 ) + A^ p 1 ^ 2 dp = 0. Equation (12.51) 
implies that C u u — -yd(|i/| 2 ) + A'yp^^dp = 0. It follows 

C u u-d{\u\ 2 ) =0. (3.10) 

Direct computation in a local coordinates shows that 

dh = \vkmdx k A dx m , 

which is by definition the vorticity (2-form) of the flow. However, 

C u u - d(|it| 2 ) = (u m d m (u l Gik) + u l GijdkU?^jdx k - (u l Gijd k u j + u j dk{GijU l )^jdx k 

= u m [d m {G ik u l ) - d k (G im u l )\ dx k 
= —u m L 0 kmdx k , 

so by (13.101) one infers (13.91) as desired. □ 

The following is the second main result of this paper. 

Theorem 3.1. The only C 2 (M) D C l (M) solution to Problem (SI) which satisfies pu° > 0 in M 
is the special subsonic solution U_. 

Proof. 1. Since E and s are constants in A4, we could use Lemma 13.21 For our special case that 
AA = (ro,ri) x S 2 , we have k = 0,1, 2 . We also assumed that u° > 0. So on Mi, from the boundary 
condition u' = 0 , (i.e. u 1 = 0 , u 2 = 0 ,) there holds 

2 

lo\2 = 5>(G^) — d2(Gnu l )\ = d\(Go2U°) — ^(Goin 0 ) = 0. 
i=0 

The last equality follows from the fact that Gij = 0 for i j. Hence we conclude that for all m,k, 
cOmk = 0 on My; namely, dfZ = 0 on M\. 

2. Now acting d on equation (13.101) . we have C u du = 0 in JA. By Lemma A.5 in [7J p.2539], 
one concludes that d-fZ = 0 in M. Since the first Betti number of M is zero (that is, M is simply 
connected), by Poincare lemma, u is exact: there is a function ip € C 3 (M) so that dp = u, or 
u = grad<^. So conservation of mass becomes 

div (/9 • grady?) = 0 , 


(3.11) 






THREE-DIMENSIONAL STEADY COMPRESSIBLE EULER SYSTEM 


11 


with p a function of |grad<^| 2 determined by the Bernoulli law ^|grad(/?| 2 + = E. Equation 

(13.111) is the so-called potential flow equation. It is a second order elliptic equation for subsonic 
flows. 

3. The boundary conditions for the potential function ip are 

ip = 0 on Mi, (3-12) 

\ gra dip\ 2 = K 0 A 2 (E-^(^)^ on M 0 . (3.13) 

The theorem is proved by applying the following Lemma 13.31 □ 

Lemma 3.3. Let ip € C 2 (M) be a solution to problem (13.111) (13.131) so that d r ip > 0 on M$. Then 
it must be a spherically symmetric subsonic flow. 

Proof. 1. We use the standard global Descartes coordinates of M. given by (z 1 , £ 2 , z 3 ) so that 
r = x° = \J ( z l ) 2 + (z 2 ) 2 + (z 3 ) 2 . Then one checks that equation (13.111) could be written in the 
non-divergence form as 

3 3 

Y a lj (Dip)dijip = c(\Dip\ 2 ) 2 Ap - Y dupdjipdijip = 0. 
i,j=1 i,j=1 

Here and in the rest of the proof, di = dij = qTq z . , and c is the sonic speed; Dip and A ip are 
respectively the standard gradient and Laplace operator in R 3 acting on ip. 

2. Let ip s = ip s (r) be a special symmetric subsonic solution to problem (13.111) (13.131) . The 
existence of ip s is guaranteed by Lemma l3Jl We see that w = (p s — ip solves the following equation: 

3 3 

a l i(x)dijW + ^ b t (x)diW 

i,j=1 »=1 

3 3 

= Y a 12 (Dip s )w + Y dijV ' ( alj ( D Ps) ~ a lj (Dip )) = 0. 
i,j= 1 i,j= 1 

Since ip s is a uniformly subsonic flow, this is a linear uniformly elliptic equation of w in A4. 

The boundary conditions of w are 

w = 0 on Mi, Dw ■ l = 0 on Mo, 

with l = Dip s + Dip. By the assumption that d r ip > 0 and the fact that d r ip s > 0, we see l ■ 
(r°, 0,0)/|r°| = d r ip s + d r ip > 0, hence on Mo we have an oblique derivative condition of w. By 
maximum principles m Theorem 3.5 in p.35] and Hopf boundary point lemma [T7, Lemma 3.4 in 
p.34], we conclude that w = 0 as desired. □ 

Remark 3.2. The uniqueness of symmetric subsonic solutions to potential flows in infinite conical 
nozzles had been proved in [22| by applying Harnack inequalities, see also m for the uniqueness 
result of the case that M = T 2 . The existence of isentropic irrotational subsonic flows in general 
three-dimensional largely-open nozzles was proved in [23j, while the same existence problem for the 
three-dimensional Euler system still remains open. See [6] and references therein for some results 
on subsonic flows in strip-like domains for the two-dimensional compressible Euler system. 
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4. Stability of spherically symmetric subsonic flows 

In this section we continue to investigate Problem (S). We are wondering if a spherically sym¬ 
metric subsonic flow Ub (called as a background solution in the sequel) constructed by Lemma 13.11 
is stable under multidimensonal perturbations of boundary conditions. The following is the third 
main theorem we will prove in this paper. 

Theorem 4.1. Suppose that Ub is a background solution so that for x° G (ro, r \) and t = M 2 (x°), 
there holds 

7(1 + 2 7 )f 4 + (—4 7 2 + 2 7 - 3)t 3 + (14 - 77 )t 2 - 19* + 6 < 0; (4.1) 

here Mb = Ub/cb is the Mach number of the background solution. Then for a G (0, 1), there exist 
positive constants £0 and C depending only on the background solution Ub and ro, 77 , 0,7 so that if 

ll U l( X ) IIC 3 ’“(Mi) “I" ll^ 1 ^ ) _ Eb \\c 3 ’ a (M 1 ) ll Sl ( X ) ~ Sb \\c 3 ’ a (M 1 ) 

+ \\po(x') - Pfe|| C 3, a(Mo) < £ < Co, (4.2) 

there is uniquely one solution U to Problem (S), with p G C 3,a (M), s,E,u G C 2,a (A4), and 

\\P ~ Pb\\c3.<*(M) + II s — s b\\ C 2 -“(AT) + \\ E ~ E b\\c 2 ^(M) + IKIIc 2 .“(AT) — ^ £ ‘ ( 4 - 3 ) 

Remark 4.1. The requirement (14.1j) is a quite rough condition for the stability of Ub derived from 
the decomposition stated in Theorem l2.ll (cf. Remark 14.21) . Note that /(1) = — 2(7 + l ) 2 < 0, so a 
sufficient condition for m is that the Mach number at the entry ro is quite close to 1 and 77 is 
close to ro- 

To prove Theorem 14. 11 firstly we formulate a nonlinear problem (S2) by specifying the functions 
L k ,Lk appeared in Theorem 1 2. 11 and by the same theorem, we infer that Problem (S2) is equivalent 
to Problem (S). Then we construct a nonlinear iteration mapping and solve Problem (S2) by using 
a Banach fixed point theorem. Lots of expressions derived here will also be used in Section [5] for 
the studies of transonic shocks. 


4.1. Problem (S2). To formulate Problem (S2), we need to compute the exact expressions of 
(12.1311 and (12.1111 and then specify the auxiliary functions L k ,Lk- 


4.1.1. Specification of boundary conditions. By the definitions of divergence operator and covariant 
derivative in a local spherical coordinates system of A4, we have 

2? ,0 1 / s o?/° 

divu = d 0 u° -\ -p- + —j=d a (VGu 01 ) = d 0 u° -|-g- + d a u a + u 1 T\ 1 , (4.4) 

X U yj G ' ' XT 

D u u + -gradp = (v?djU k + u 3 u m T k rn + - dipG lk )dk■ (4.5) 

P P 

It follows, by setting Ls(-) = — 9q), the identity 


,D u p 


grad p 


TP P 

+ ^ - %G a pu a u^ + 


— 1 0Q P + UU 


pu °(—— + div u + Ls(D u u + L)) = 

“ \ c~ j 

ApE 2p„ n R 1 u a “■ 5 


(u°) 


0\2 


0..Q; 


1 1 

+ 




n 47 1 

OaP - T-rjP 

7 - 1 r 


d a A(s)p^ + ^d 5 (G a pu a uP) - ^ d a E 

+ ^d a (VGu a ) - pu a uPr%. (4.6) 


— 1 u° 
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Comparing this to (12.1311 , we see it is actually a nonlinear Robin condition for p on M \: 

1 2y p(u 0 ^ 2 


doP + 


x° ( u °) 2 — c 2 


= G 1 (U) + G 2 (U), 


with 


Gi 4 


G 2 4 


+ 0 ) 2 i 

“72- 1 


-^d a (VGu«) 


(4.7) 


(4.8) 


1 u c 


(«o )a _ 1 L 7 _ ! - ^(G^«V) + 


ir 


+pu < Vr° /3 -u%° d aP 


(4.9) 


4.1.2. Specify equation of pressure. We now compute the explicit expression of the equation (12.111) . 
It is straightforward to check that 


D u (M). div ( _ clcl(Du ® Du) = 2- 
V 7P / V P J " IP 


+ 


1 

7P 


u 0 ) 2 - c 2 


((u 0 ) 2 - c 2 ) d 0 p + u°d 0 u°d 0 p - ^-(dop) 2 + ^d 0 pd 0 p 


)doP 


(x°y 


r A'p 


P 


P 


u 


0 \ 2 


-(<9 0 u u ) 2 + 2 0 

where A' is the Laplacian on S 2 , and 


+ 2—<9 0 u° + Hi(U), (4.10) 


tfl(tf) = —-“ + 2 

x u p z 


11^ . 1 - 

- djkP H- d k u J djp - g u k u J dkpdjp 

7p 77 ? 7 


A_ 2 7( — +divtt| . (4.11) 


(fcj)#(0,0) 

— ^ ( dju l diu? + 2 Tjpu^diu^ + 

0 j)tRo,o) 

Replacing terms like doit 0 , (<9ou 0 ) 2 , (u 0 ) 2 in (|4.10l) by using suitable Euler equations, after some 
straightforward computations, we get the identity 


x u \ 7 p 


d„ i Ah') - div (+2lA -c;7(du®d») 


1 

7P 


7P 


2E- 


7P 


\ P J 


7+1 


-c 2 )^p 




(x 0 ) 2 
■ 4 1 


A'p 


+ 


7pa: u 


£?- 


7 c 

7-1 


<9oP 


7-1 4 7 £| _ _£L 

7-1 / 

+#i + H 2 + R 3 , 


(<9op) 2 + 


(x 0 ) 2 


E — 


7-1 


(4.12) 


H 2 = — 

IP 


d 0 pG(D u u + (grad p)/p, d Q ) - p 7 l d 0 p-^D u A(s) 


-d 0 p(u a d a u° + u a uPr° a/3 ) + p 1 1 d 0 p^d a A(s) 


with 
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( ~ 7 T ) G(D U U + (grad p)/p, 0 O ) 2 + 2d °^ G(D u u + (grad p)/p, 0 O ) 


p{u°) 


+ -j-^{u a d a u° + u a u^T° a p)G{D u u + (grad p)/p, 0 O ) 

(u ) 

- ( u a d a u° + u a u^T° a p) 2 - -^^(u a d a u° + u a uPr° a/3 ) 

+—qG(D u u + (grad p)/p,d 0 ) - —(u a d Q u° + u a u p T° a p), 

X u Xq h 

r i 2 2 

H 3 = -G a pu a u 13 —0q P~\ - Q d 0P + Tl)Y 2 

I'lP 7 px u (x v ) z 


(4.13) 


+ 


(do pf 
7 p 2 


-1 + 


c 4 1 


1 2E-^ l 2E- G a pu«uP - ^ 


n 2c 2 

7 “ ' 


(4.14) 


Now set 


o ?/ o / xy-i \ 

L °(^) = L\<p 3 ) = 0, T 2 (+ 4 ) = ++ 

£ 3 (+p) = - f— + + T-7)Y2( ua da u0 + + 4, ) G((p 0 ,d 0 


+ (^oy2 G, ^°’ 5 °) 


7P ' p(u 0 ) 2 (u 0 ) 2 

2 


By multiplying 7 p to both sides of (I4.12p . and comparing it to (12.111) . we see (12.111) is equivalent 
to the following second order equation of pressure 


N(U) = 


7 + 1 

-( 

7-1 


2 E - ---c“ 

,2 „4 


! )+> 


(x°y 


rA 'p 


+ ^\E- 


— \E- 
P V 


+ 


x u 

\2 


7 c 

7- 1 


0 O P 


7-1 4 7 E - -7 

7 — 1 / 


^0 P) 2 + ^ 




7-1 


= +)([/), (4.15) 


where 

— F\(U) = 7pi?3 + 


-d 0 p{u a d a u° + u a vPT° a p) + p 1 1 d 0 p^d a A(s 


-IP 


(2 q) ( u a d a u° + u a u p T° ap ) + 


20 o p 2 


+ — 


(■ u a d a u° + u a vP r° a p) 


p(u 0 ) 2 Xq 

+u a d a p + ^G a/3 d a pdpp + ^ f u k u j djkP + u k dku j djp - -u k u j dkpdjp 

X P , , n nl V P 


(fc,j)+(0,0) 

-7P ^ ( dju l diu j + 2 T l j/ 3 u l 3 diu j + 

Gj)+(o,o) 


(4.16) 


4.1.3. Problem (S2). From the above computations, by Theorem 1 2. 11 we see that Problem (S) could 
be written equivalently as the following Problem (S2), for those unknowns p,p,u with regularity 
as assumed in Theorem 12+1 

Problem (S2): Solve (12.91)(12. 101) (12. 121) and (14.151) in M, subjected to the boundary con¬ 
ditions (13.11) (13.21) and (14.7(1 . 
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4.2. Problem (S3). Since we are dealing with a small perturbation problem, in this subsection 
we separate the linear main terms from the nonlinear equations (14.71) and (|4.15|) , and write them in 
the form C(U — Ub) = M{U — Ub), where C is a linear operator, and N(U) are certain higher-order 
terms defined below. By this way we formulate Problem (S3), which is equivalent to Problem (S2). 

Definition 4.1. For U = p,E,p,s,u etc., set U = U — Ub . A higher-order term is an expression 
containing either 

(i) U\dM and/or its first-order tangential derivatives; 
or 

(ii) product of U, and/or their derivatives DU ,D 2 U , where D k u is a k th -order derivative of U. 


4.2.1. Linearization of boundary conditions. Recall that for the background solution, pb solves (|3.5|) . 
so (14.7p is equivalent to 


do (p~Pb) + 


27 


p(u 0 ) 2 
(u 0 ) 2 — c 2 


Pb{u° b ) 2 \ 

K) 2 -c 2 J 


g 1 + g 2 . 


(4.17) 


Using expressions like 


c 2 - cl = 1 - ~(p~Pb) + Pb 1 {M S ) ~ A ( s b)) + 0(l){p - p b ) 2 + 0(l)(A(s) - A(s b )) 2 , (4.18) 

Pb 

where 0(1) represents a bounded quantity with the bound depending only on the background 
solution Ub and \U — £/&|, after some straightforward computations, (14.171) could be further written 
as 


doP + 7iP = G = Gi + G 2 + G 3 , 

with 7 i a constant determined by the background solution at x° = r\\ 

2 7 M fe 4 -M 2 + 2 


and 


G, = - 


2 7 


K) 


0\2 


A 

7i = 


ut 


x° ( M 2 - l) 2 > °’ 

u 2 (c 2 - c 2 ) - c 2 ((u °) 2 - u 2 ) 


(■ u°) 2 — c 2 u? — ci 


2 -72 P + Pb 


'll 2 — c 2 
u b c b 


(4.19) 

(4.20) 

1 


2 b ° b 2 ,M E ~ E b) ~ G a0 u a uP] + 

\ u b c b> 


2 Pb E b 

(“? - 4Y 


+ pr\A{3)-A( Sb )j \\. 


(■ u°) 2 — c 2 u 2 — c 2 

G(l)(|p | 2 + |A(s)| 2 ) 

(4.21) 


We note that G 2 , G 3 are higher-order terms (the terms with underlines are given by boundary data, 
so fulfill the item (i) in Definition 14.11) . while Gi depends on the boundary values of u 1 ,u 2 on M\. 
Boundary condition (I4.19P is an equivalent form of (|2.13D . 


4.2.2. Linearization of pressure’s equation. For (14.151) . note that N{Ub) = 0, we may get a linear 
operator L and write L{U) — F 2 (U ) = N(U) — N(Ub), with F 2 (U) a higher-order term. Then (14.151) 
becomes L{U) = F\{U) + F 2 (U). In fact, using (14.181) . by setting t = u 2 /cl = M 2 € (0,1), direct 
computation yields that (14.151) can be written as 


L{P ) - --rov2 A 'P + ( t ~ l ) d oP+ -^K^doPF TU)Y2 e ( t )P + r7F2 rf iW^ + 


[x 


m 2 




(x‘ 


0^2 


(X' 


0)2 


7 

Pb 

0)2 


(x 1 


d 2 (t)A(s) 
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with 


m = 


(4.23) 

e{t) E 

(t _ 1)3 N 1 + 27 ^ 4 + ( 4 / 2 + 2 7 3 )t 3 + ( 14 7 t)* 2 19t + 6 ] > 

(4.24) 

di(t) 4 

(t _ 1)3 ((27 3)t 2 + 8t 3), 

(4.25) 

d 2 {t) ± 

7 /i ( t-i)3 [2 +(y ~ - 3 >‘ 2 + 8t - 3 i. 

(4.26) 


and 

~F 2 = 


4:^P 

( x °) 2 


E - 


1 


4 - - / A 7 / 2 2 


1 


"( C - C b)pj2 A '^ + 


2 i? — 


( E - 


(c^-c) + „ 

7—1 / x u \ 7—1 


(c 2 - ci) 


—-{do pf + d 0 {p + Pb)doP 
Pb 


7- 1 

Pb p 


dip 


{dp Pbf 
Wb 
{doPbf 


<? + <% 

P 

1 1 




|tt| 

p 


ur 




Pb 


Pb \p Pb 


!PbU% 


<{dppf 


+ 


7 V l M l 
47 p b 


ut 


-M + 


2 c? 


P PbJ PbU z b V 7 - 1 


c“ - cr 


— E 


do(p + Pb)d Q p ( _c/_ 

7 VpM : 




47 9 0 p b 7 + 1^2 , 2 (5oP&) 2 

-9 0 Pfe + 


7 - 1 Pb 


(i + E-5 


1PbU b 


7 — 1 (x°y 7 — 1 x° 7 — 1 
x [ 0 (l)p 2 + 0(1)(j4(s) - A(sfe)) 2 ] . 

We easily see that (J4.22I) is an elliptic equation of (perturbed) pressure. 

4.2.3. Problem (S3). We now reformulate Problem (S) equivalently as the following Problem (S3). 


7-1 u 2 b/ 

(4.27) 


Problem (S3): Solve functions U = (E, A(s),p,u' = u a d Q ) that satisfying the following 
problems (14.28[> — (14.31 j) . 


D U E = 0 
E = E l {x') 1 

D u A{s ) = 0 
/4(s) = ^(si)( 

' L{p) = F(U) 


' 1 ? 

Mr, 

(4.28) 

in A4. 

- i 

on Mi ; 

(4.29) 

in M., 

on Mo, 

(4.30) 

on Mi; 

,) = 0 in A4, a = 1,2, 

on Mi. 

(4.31) 


It/ = u\ (x ! ) 

We note that (I4.30P is a mixed boundary value problem of a second order elliptic equation, while 
(14.281) . (14.291) and (14.311) are Cauchy problems of transport equations. In fact, we could write the 
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equation in (14.311) in a local spherical coordinates as 


1 


uhd jU a + u^u m T^ m = —d iP G ia . 


P 


For a = 1, it reads 


with 


27 /° 

v?djU 1 H-u 1 = — 


(x°) 2 p 


dip + hi(U), 


hi(U) = (sinx 1 cosx 1 )(u 2 ) 2 H—— rt°)u 1 . 


(4.32) 

(4.33) 

(4.34) 


For a = 2, it reads, in the local coordinates used above, that 


2 1 
u^djU 2 H— ttU°u 2 = —2cotx 1 u 1 u 2 — 


1 


yd 2 p- 


x” / o(x°sinx 1 ) 2 

This equation has an artificial singularity when sinx 1 = 0. It is obvious that we can avoid this by 
using another local (spherical) coordinates. So by symmetry of u 1 and u 2 , it suffices to solve (I4.33P 
and obtain an estimate. 

4.3. Proof of Theorem 14.11 For k = 2,3, and U = ( p,s,E,u' = u l d\ + u 2 d 2 ), suppose that 
p G C k,a (M ) and s,E,u^dg € C k ^ 1,a (M). We define the norm 

2 


k — \\p\\c k ’ a {M) + ll' s llc" c - 1 ’“(A!) + \\^\\c k - 1 ’ a (M) + ^2 

13=1 


fi 


C fc -L“( M) 


(4.35) 


Let K be a positive number to be chosen. We define Xk e to be the (non-empty) set of functions U 
so that (13.11) and (13.21) hold, and \\U — £4|| 3 < Ke. To prove Theorem 14.11 we construct a mapping 
T on Xke for suitably chosen K and e, and show that it contracts under the norm Then by a 
Banach fixed point theorem, T has uniquely one fixed point U G X^e, which is exactly a solution 
to Problem (S3). 

4.3.1. Construction of mapping T■ For any U G Xke> by the following process we define a mapping 
T : U i-A- U. Set U = U — Ub. Then we only need to determine U. We also use C to denote generic 
positive constants which might be different in different lines. 

Determination of E and A(s). Noting that $i, and Eh are constants, we solve the unknowns 
E, A(s) from the following Cauchy problems of linear transport equations, where the velocity field 
u is given as part of U G Xk £ '- 


D U E = 0 in At, 
E = E\ — Ef, on M \; 


D u A(s ) = 0 in M, 

A(s) = A(si) — A(sb) on M\. 


(4.36) 


Since u G C 2 ’“(A4), and recall that E\ — £),, si — Sb G C 3 , “(Mi), by Theorem 16.11 we have 
Lemma 4.1. There are uniquely E,A(s) G C 2,a (Xi) solve (I4.36p . In addition, 


E 

+ 

A& 


C 2 ’ a (M) 



_ < Ci 
C 2 ’ a (M) 


E 


C 2 ’ a {M i) 


+ 


Ms) 


C 2 ’ a (M i) 


<C- l£ , 


(4.37) 


with C\ a positive constant depending only on Ub- The second inequality holds provided that (|4.2I) 
is valid. 


Once we solved E and T(s), we get E = Eb + E, A(s) = A(sft) + A(s). 
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Determination of p. Now consider the following mixed boundary value problem of p (comparing 
to problem (14.301) 1: 

£(p) - - + (* - l)fio# + + j^yP(t)p 

= + nv\ (4 38) 

p = p 0 -p b on M 0 , 
doP + 71 P = G(U) on M\. 

Note here that for the two terms d\{t)E and ^S^ d2(t)A(s) on the right-hand side, we take E 

and A(s) to be the functions solved from Lemma 14.11 

We now specify the nonhomogeneous term F{U). In F-i (see (14.271) 1. all U should be replaced 
by U — Ub (recall that the U € Xx e has been fixed). So by direct computations we get that 
11 F 2 11 c i , a < CK 2 e 2 . In the expression of F\ (see (14.161) 1. we take all U to be the given one. So 
by the smallness of u a , we have < CK 2 £ 2 ,\\Fi\\ cl , a Q^^ < CK 2 e 2 . Therefore we 

obtain, for a positive constant C depending only on Ub, that 

\\F(U)\\ c i, a(M )<CK 2 e 2 (4.39) 

Next we specify the term G(U ) in boundary conditions. For G\ (see (14.8p ). the u a should be the 
boundary conditions (u\) a (hence belong to C' 3 , “(Mi)), and the others are replaced by the given 
U. So by (14.21) we have || (^1IIc' 2 ,“(a^ 1 ) — Ce. For the term G2 (see (14.91) 1. u a and A(s),E should 
be the given boundary data (u\) a and yl(si),i?i respectively, while the others are replaced by the 
given U € Xk £ ■ Hence it still lies in C 2 ,ol {M\ ) and we have ||G f 2|| C 2,a(jvf 1 ) < Ce 2 + CKe 2 . For the 
expression of G 3 (see (|4. 211) 1. except the underline terms are replaced by the given boundary data, 
all the other U can be taken as the given U in Xke, and it easily follows that \\G^,\\ C 2, a ^ M ^ < 
CK 2 e 2 + C£. So in all, we obtain 

\\G(U)\\ c2 , a{Mi) <CK 2 £ 2 + Ce, (4.40) 

if we choose later, without loss of generality, that K > 1. 

Under the assumptions of the Theorem 14.11 the coefficient e{t) is nonnegative. Recall also that 
71 > 0 (see (14.201) 1. So by the standard theory of second order elliptic equations with Dirichlet and 
oblique derivative conditions (Em Chapter 6 ]), we have 


Lemma 4.2. There is uniquely one solution p £ C 3,Q (AI) to problem (I4.38j) . In addition, there is 
a constant C depending only on the background solution Ub and Xi, so that there holds 


II C 3 ’ a (M) — 






+ 






+ 


E 


< C(K 2 £ 2 + e). 




(4.41) 


Hence we obtain that p = p + pb- 


Remark 4.2. We see that (14.11) is used to guarantee that problem (|4.38[) has uniquely one solution. 
A detailed study of the spectrum of the operator £ with the homogeneous boundary conditions 
(like what we do in Section 5 by separation of variables) would definitely refine it. However, for 
our present purpose of clarifying the basic ideas, we are content ourselves with (|4.1D . 
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Determination of u'. Now we solve u a (a = 1,2) from the transport equations G(D u u + 
gl ^ d:p , d a ) = 0. For a = 1, we have the problem 


u j djU 1 + = -jjpjtjdi p + hi(U) in M, 

u l = K) 1 on Mi. 


(4.42) 


Here [u\) 1 is the given boundary data; p on the right-hand side is given by Lemma 14.21 while u 
and U in 

hi(U) = (sinx 1 cosx 1 )(n 2 ) 2 + (u b — u°)u 1 (4.43) 

are the one we had fixed in X^e- 

We have the following lemma due to Theorem 16.11 


Lemma 4.3. There is uniquely one solution u 1 € C 2,a (JA) to problem (14.421) . In addition, 

ll^ 1 llc 2 ’“(Al) — ll U l 1 |lc 2 .“(Mi) + C ( Mc^(M) + \\M U )\\c^(M) ) < c(k 2 E 2 + E^J (4.44) 
for a positive constant C depending only on the background solution and ro,ri. 


As explained earlier, we could change the coordinates to solve u 2 and similar results hold. Then 
from U € Xk e we obtain uniquely one U = E,u a d a ), and the estimates (14.371)(14.411)(14.441) 

< C(K 2 e 2 + e). By choosing K = max{2C, 1}, and £q < min{l/iL 2 ,1}, 


yield that 

we have 
T:X Ke 


U-U b 

u-u b 


< K£ for all £ < £q. Hence we proved that U € X^ e , and the mapping 


Xke is well-defined. 


4.3.2. Contraction of the mapping T. For U^ l \U^ € Xk £ , set U^ k] = T(U^), k 
to show that if £q is further small, then 


t/ (1) - f/ (2) 


< - 

2 “ 2 


[/(L _ [/(2) 


1,2. We wish 


(4.45) 


The idea to establish (14.451) is to consider the problems satisfied by U = — U^ 2 \ Since the 

process is standard once we understand the definition of T, and is quite similar (but simpler) to 
that described in Section 15.8.31 we omit the details here. 

Then by Banach fixed point theorem, T has one and only one fixed point, say, U, in Xk £ ■ By the 
construction of the mapping T, the fixed point is a solution to Problem (S3). On the contrary, for 
a solution to Problem (S3) which lies in Xx e , it must be a fixed point of T. The proof of Theorem 
l4~H is completed. 


5. Stability of spherically symmetric transonic shocks 

In this section we study the stability of spherically symmetric transonic shocks under multidimen¬ 
sional perturbations of boundary conditions, which was treated in jl] by using the “non-isentropic 
potential flow model”, and in [7] by considering the full Euler system. In [7] only uniqueness was 
proved; namely, if the perturbations of the upcoming supersonic flow and the back pressure are 
small, then there will be only one transonic shock solution in some function space, provided that 
it exits. We will show below that in a suitable class of functions, a disturbed transonic shock solu¬ 
tion does exist and is unique in a neighborhood of the background solution U b , provided that the 
S-condition, which also occurred in [7], is valid. The main difficulty is to decompose the Rankine- 
Hugoniot conditions (R-H conditions) on shock-front in a way compatible to the decomposition of 
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the Euler system stated in Theorem 12.11 We will discover many subtle intrinsic structures lying 
beneath the related free boundary problem. 

In the following we firstly review the formulation of the transonic shock problem (T) and the 
existence of background solutions, and state the main theorem, i.e. Theorem 15.11 Then we start 
to decompose the R-H conditions, and formulate Problem (T) step by step into the more tractable 
problems (Tl), (T2), (T3), (T4), each of which is equivalent to Problem (T). Finally, Theorem 15.II 
is proved by applying a nonlinear iteration method to the nonlinear free boundary problem (T4). 

5.1. Transonic shock problem (T) and main result. Let Ml = (ro,ri) x S 2 be a spherical 
shell with entry Mq = {ro} x S 2 and exit M\ = {ri} x S 2 . We use U = (u,p,s) to represent the 
state of the gas flows in Ml. Suppose that 

= {(x°, x') €M:x°= ^(x'), x' G S 2 } (5.1) 

is a surface, where ^ : S 2 —>■ Ml is a C 1 function. The normal vector field on is given by 

n = (d a ^G^dp-G 00 d 0 )\ s (5.2) 

We also set = {x G Ml : x° < ip(x r ), x' € S' 2 } to be the supersonic region , and M^ = {x € 
Ml : x° > ip(x'), x' € S 2 } to be the subsonic region. 

Definition 5.1 (Transonic shock). Let ^ G C 1 (S 2 ) and U ± G C l (M.±) n C(M^). We say that 
U = (U~, U + ;ijj) is a transonic shock solution , if 

1) L/ ± solve the Euler system (11,41) in M.± in the pointwise sense; 

2) U~ is supersonic, and JJ + is subsonic; 

3) the following R-H conditions hold across S^\ 

[[G(u, n)pu + pn ]] = 0, [[G(u, n)p ]] = 0, [[G(u, n)pE]\ = 0, (5.3) 

where [[/(H)]] = f(U + \ s ^) — ) denotes the jump of a quantity f{U ) across S^; 

4) there holds the following physical entropy condition 

[[?]] =P + \s^ ~P~\s* > °- (5- 4 ) 

By the definition we infer that a transonic shock solution is a weak entropy solution of the steady 
Euler system (c/. Section 4.3 and Section 4.5 in [13] ). 

To formulate the transonic shock problem, we also need to specify boundary conditions. Since 
the flow U~ is supersonic near the entry Mq , the following Cauchy data should be given: 

U = Uq {x') on M 0 . (5.5) 

Here we also require that (if°)o > c o t° m &ke sure the steady Euler system is symmetric hyperbolic 
in the positive x°-direction on Mq. On the exit Mi, it is known now that one and only one boundary 
condition is necessary. From the physical considerations, as in the studies of subsonic flows, we 
propose the pressure 

p = pi(x') on Mi. (5.6) 

It turns out that this is also a mathematically challenging boundary condition ( cf. [30] and discus¬ 
sions at the end of [29]). 

Problem (T): Find a transonic shock solution in Ml which satisfies the boundary condi¬ 
tions (15.5|) and (15.61) pointwisely. 
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5.1.1. Background solution. The existence of spherically symmetric transonic shock solutions U b = 
(Ufl, Uj(~; 77 ) (which are called as background solutions in the sequel) to Problem (T) has been 
established in m (see also m)- For given [ro,ri], we note that Uj^ is determined by, and actually 
depends analytically on, the five parameters ( 7 , 77 , 7 ^( 77 ), 75 ^( 77 ), £( 77 )) with t(r b ) = (M(f) 2 (rb) € 
( 0 , 1 ), 7 > 1 , 77 € ( r °, r 1 ) j p(f (r b ) > 0, and p^(r b ) >0. As a consequence, all the coefficients 
on the left-hand side of (I4.22P . such as e(t), pflhfl), depend analytically on these parameters. In 
addition, it is useful to note that the spherically symmetric subsonic flow U(f could be extended to 
[77 — hf,r{\ x S 2 , for a small constant hf depending solely on these parameters ( cf. (31,, p.323]). 


5.1.2. The main result. We now state the forth main result of this paper. 


Theorem 5.1. Let U b and 7 satisfy the S-Condition (see Definition 15.71 and Lemma \5.f\ ), and 
a € (0,1). There exist £0 and C * depending only on Ub and 7, a,ro,ri such that, if the upcoming 
supersonic flow U(f on Mq and the back pressure on Mi satisfy 

IK - U h Hc 4( Mo ) — £ — £ 0; (5.7) 

\\pi-Pb\\c3, a{Ml ) < £ < £ 0, (5.8) 

then there exists a transonic shock solution U = (U ~, U + ] if) to Problem (T), so that if € C' 4 ’“(£ 2 ), 
U- € C 4 (MT), p + e C 3 ’“(AI+) ; u + ,p + ,s + € C 2 ’ q (M+), u + \ s ^p + \ s ^s + \ s ^ € C 3 ’“(S 2 ), and 



IIV 7 — r fellc 4 ’“(s 2 ) — C* e i 

(5.9) 


W U ~ U b \\c*(Mj) - C * £ ' 

(5.10) 


II 77+1 — 77+1 II 4- II/7+ _ 77+11 < n <r 

II u Is7 u b l5V'llc3.-(S 2 ) ^ H U U b II 3 — Ly * fc " 

(5.11) 

Furthermore, 

such solution is unique in the class of functions if,U~,U + with 



IIV 7 - ^ll C 3.-(5 2 ) ^ c '* e > 

(5.12) 



(5.13) 


II 77+1 — 77+1 II 4- 1177+ _ 77 +II < n r 

II u Is7 u b lsV> Ilc 2 , a ( 52 ) + II u u b II 2 — Ly * fc " 

(5.14) 

The norm |- 

1/. is defined bu (14.3511. with A4 there revlaced bu A4 + . 



5.1.3. Existence of supersonic flow. The existence and uniqueness of supersonic flow U~ in M. = 
(ro,ri) x S 2 subjected to the initial data U(f satisfying (15.71) follow from the theory of classical 
solutions of the Cauchy problem of quasi-linear symmetric hyperbolic systems if £q is sufficiently 
small and |ri — 77 1 is not large (cf. J2j). Furthermore, one can obtain that 

II ^ — U b ||c'4(3vt) — ^ l£ ’ (5.15) 

where C\ > 0 and eo > 0 depend solely on Uff (rfl) and 77 — tq. This implies (15. lOf) . So Problem (T) 
is actually a one-phase free boundary problem, for which the free boundary (i.e. the shock-front) 
and the subsonic flow U + are to be solved. For simplicity, from now on we write U + as U. 


5.2. Reformulation of R-H conditions. For transonic shock problem, the R-H conditions repre¬ 
sent a nontrivial nonlinear structure of a discontinuous flow held. In this subsection we decompose 
them to find their roles in determination of shock-front and subsonic flows. For the sake of com¬ 
pleteness, we shall take some computations from [7J pp.2520-2523]. 
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5.2.1. Decomposition of R-H conditions. Consulting (15.21) . we see the normal 1-form on the shock- 
front (15.11) is given by h = dip — da; 0 , and the R-H conditions (15.31) could be written in terms of 
differential forms as 


\[n{u)pu + pn]\ = 0, 

(5.16) 

[[n(u)p]] = 0, 

(5.17) 

[[n{u)pE}\ = 0. 

(5.18) 


Here n{u ) is the action of the 1-form n on the vector field u. 

By the definition of shock-front, the mass flux m = — n(u)p \ s </, = — n(u~)p~\ s ^ 0. So from 

(15.171) and (15.18|) we infer that E~ \ s ^ = E + \ s ^. Combining this with (12. 1|) , one concludes that the 
Bernoulli constant E is invariant along flow trajectories, even across a shock-front; so it is actually 
determined by the upcoming supersonic flow Uff. 

Set u = uq + u' = u° do+ u a d Q , u = u° + u' = u°Gqo d x° + u a G a p dx@. Then to (it) = d ip(u') — u°, 
and m = pu 0 — d ip(pu'). So ()5. 171) and ()5. 181) may be written as 


while (15.161) is decomposed to be 


[H] = o, [[E}] = 0, 

[[mu°+p\] = 0, 
[[pdip — mu']] = 0. 


(5.19) 

(5.20) 

(5.21) 


If [[p]] > 0 (which is guaranteed by the physical entropy condition satisfied by the background 
solution, and the small perturbation estimate (15.111) to be established), from (15.211) we solve that 


A m[[u']] 
dip = u = 


M] 


Si’ 


= Rofi*(u') + g 0 , 


with fiQ a positive constant, and go a 1-form on S 2 : 

(pu°)i 


Ro = 


Pb ~Pb 


> 0 , 


n 


go = go(U,U~,Dip)± 


rrofftt 7 ]] 


M 


Si’ 




(5.22) 

(5.23) 

(5.24) 


Here we use the pull back if*, so if*(u') = u'\sii is a 1-form on S 2 . We note that go is a higher-order 
term (see Definition 15.21 below). which depends on U\si>, and Dip. 

The R-H conditions (I5.16l) - (l5.18p are equivalent to (15.191) (15.201) and (I5.22|) . if [[p]] / 0. 


Remark 5.1. By (15.221) . it is necessary that doj = 0, which implies 

d(ip*(u')) = — -dg 0 . (5.25) 

Po 

On the contrary, since the first Betti number of S 2 is zero, by Poincare lemma, (|5.25D is also 
sufficient for the existence of a function if v on S 2 so that (15.221) holds, and f s2 ip p vol = 0 (here vol 
is the volume 2-form on S 2 ). Such a function ip p is called the profile of the shock-front ip. We also 
call the number r p = ip — ip p as the position of the shock-front. 
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5.2.2. Linearization of R-H conditions. We recall the following definition from [7J p.2522]. 

Definition 5.2. Set U = JJ + — U^~ and f = r p — r b . A higher order term is an expression containing 
either 

(i) U~ — Uf and its first or second order derivatives; 
or 

(ii) the products of ip p , r, U, and their derivatives DU, D 2 U, Dip, D 2 ip, and D 3 ip, where D k u are 
the fc-th order derivatives of u. 


As shown in 0 p.2523], we can write (15.191)(15.2011 equivalently as 

^*(«°) = hi iV + r p - r b ) + gi(U, U~,ip, Dip), (5.26) 

ip*(p) = 92 {V + r p - r b ) + g 2 (U, U~,ip,Dip), (5.27) 

^*{p) = 93 iV + r p -r h ) + g 3 (U,U~,ip,Dip). (5.28) 

Using A(s) = pp _1 , it follows that 

iP*(A[s)) = p A (iP P + r p - r b ) + g A {U,U~,if,DiP). (5.29) 

Here pj (j = 1.2,3,4) are constants and gk (k = 1, 2,3,4) are higher order terms. From (I5.27|) and 
(15.291) we also have 

r(A&) = -rip) +9i - -92. (5.30) 

92 92 

We then obtain the following lemma. 


Lemma 5.1. Suppose that ip and U = U + are C 1 , and the physical entropy condition holds. Then 
the R-H conditions (1 5.16 f) — (1 5.18 [) are equivalent to [[A]] \s. tjl = 0. (15.271) . (15.291) . and (15.221) . 

5.3. Problem (Tl). Now we consider the boundary condition (J2.13I) with M\ there replaced by 
the shock-front S^ and then formulate Problem (Tl), which is easily seen to be equivalent to 
Problem (T). 


5.3.1. Divergence of tangential velocity field on shock-front. By the process of deriving (14.71) and 
(14.191) . we infer the following identity 

P ™2 - + divu + L 3 (D u u + gradp ) N ) = d 0 p + 71 p - G{U), (5.31) 

/>r _ x V 7 P P J 

which holds actually at any point in A4. Now we restrict it on . So particularly x° should be 
replaced by ip. Using the commutator relation 

d a {r /) = (d a rr(dof)+r(daf), « = 1 , 2 , 


we have 


rGi = r 


= ip* 


-pu 


(u°) 


T-i 


--d a (Vgru a ) - {d a ip)r (dou a ) 


-pu 


G 0 ) 2 , 

“72 - 1 , 


(—d *(ip*u') - d a ipr(d 0 u a )) . 


Here d* is the co-differential operator on S 2 , and the last equality holds because d*u = 
a vector field u. 


(5.32) 
—divtt for 
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From (14.51) . we have 


ip*{G(D u u + (grad p)/p,G a ^dpipd a )} = {G ka (u J djU k + u 3 u m T k jm + -dipG lk )G aP dpip} 


= i/j* < 


(rt° - u a d a 'i/>)dpi/jdouP + u J u m T? m dpip + - ^ G^dppdpip > + ip* {u a )dpipd a (pp* u p ). 


s* 


4=1 




One solves ip*(dpipdouP), and (15.321) becomes 

,o \ 


ip*Gi = i>* 


pu 


(u 0 ) 2 , 

~72 - 1 , 


d *{ip*u')+ip* 


pu 


(u°) 


3 YL _ i u° — u a d a ip 


xr{G(D u n + (gradp)/p,G a %rj>d a ) - u^u m T p jm dpiP 
1 2 

— ^2 GPPdppdpip - u"^(V’V)}. 


4=1 

Similarly, one my replace the normal derivatives by tangential derivatives to compute ip*G 2 and 
ip*Gs. Now set 


T , , a 1 u 5 d 5 ip 
Li(w) = -T-fv^-n--jvr-r^ w, 


L 2 (w) = 


(u 0 ) 2 u° — u s dsip 
> 7 1 u s dsip 


(7 — l)(u 0 ) 2 u° — u s dsip 


ip*w, 


L' 3 (w) a L 3 (w) + ^ r -^jg— a(w,a ai3 9 flk d a ) + ^-jA^ G(^G s °g a gu?8s). 


Then ( 15.311) becomes 

ip* < 


—77 - f£uP + ^ivu + z}(p u E) + L 2 (D u A(s )) + L' 3 (D u u + giadp ) N ) 

«! _ ! V 7 V P ) 

= ip* {d 0 (p) + 72 P- 73 d*(^V)} + G. 

Here 72 > 0 and 73 < 0 are constants determined by the background solution, and 


G = ip* I - 


+ip* 


pu 


+ 73 d *{ip*u') + ^ip* 


^ -1 


7 1 2 — t + 2 \ 2 7 1 2 — t + 2 


ip \ (1 — t) 2 J 77 (1 — t ) 2 


ip*p 


pu 


1 


£^ 2 _ ^ u° — U a d a 1p 


+1p* 


_ 1 u ° - u a d a ip V7 — 1 


u 3 u m TP m dpip + - ^2 G^^dppdpip + u a d Q (ip*u p )dpip 


-u a d a (ip* A(s)) — pu a d a (ip* E) 
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+1p* < 


-V>* < 


+1p* < 


1 pV’ 2 ( Sr, a 0 , 2 u°g a 0U^u 5 d s (ip*u a ) 

U ds9 ^ U U + u° — u a d a ip - 


2 u° 


— I -1 


- 1 


pu a u p r° a/3 - u° 


i i 

~n + 


(u 0 ) 2 


u°dsp 


IT 


1 pip 2 0 i mra , 1 an ^ 

u -o- 9 a pu p u 3 u T m + —u d Q p 

/V\ 1 iD - u a d a ip \ K J pip z J 


+0(l)iP*{\U\ 2 + \^\ 2 + E) + -ip 


t + 


7-1 77 


(1 - i) 2 


P4 


Pfc 54 - —92 
P2 


We see () 2.13 [) is equivalent to 

d *{ip*u') = /i 5 i’*{d 0 p) + 90 V + Pe {r p - r b ) + 9 b{U, U~,ip, DU, Dp’), 


(5.33) 


(5.34) 


if we replace ip*(p) by ip via ()5.27l) . Here p b < 0, /ig > 0 are constants determined by the background 
solution, and 


1 


72 


A i 

P5 — —G H-52- 

73 73 


(5.35) 


Remark 5.2. In the expression of g b , there appear first order derivatives of p, and only first order 
tangential derivatives of A(s),u',E,ip on 5"7 Note that ()5.34l) is a first-order boundary condition 
on the shock-front. Together with (15.25j) . we have a div-curl system of the tangential velocity ip*u' 
on S 2 . 

5.3.2. Problem (Tl). For functions U = (E, A(s),p,u' = u l d\+u 2 d 2 ) and ip (note that U = U—U^ 
and ip = ip p + r p ), we formulate the following problems: 


D U E = 0 in M 
E = E~ on 5 ^; 

' L(p) = F(U, DU, D 2 p) 

P = Pi~pt 

7 * (p) = 92 {V + r p - r b ) + g 2 {U, U~, ip, Dip ); 

< 

D u A(s ) = 0 

; ip*{A(S )) = 94 {V + r p - r b ) + g A {U, U~ ,ip, Dip)-, 
dip = p 0 ip* (• u ') + g 0 (U,U~, Dip), 

d*(ip*u') = 90 ip*(dop) + 96 ip p + 96 (r p - r b ) + g 5 (U, U~,ip, DU, Dip) 
G(D u u + s y £ ,d a ) = 0 in M^, a = 1,2, 

«V = M o- 


in 7W+, 
on Mi, 


in M+, 


(5.36) 


(5.37) 


(5.38) 


on S 2 ; (5.39) 


(5.40) 


The initial data u' 0 in (15.401) is the vector field on S 2 associated with the 1-form ip*u' in (|5.39l) . 
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By Theorem 12.11 Lemma 15.II and the analysis above, it is obvious that for given supersonic flow 
U~ , the solution (U, ip) to these problems also solves the Euler system and the R-H conditions hold 
across . Hence we could rewrite Problem (T) equivalently as the following Problem (Tl). 

Problem (Tl): Find ip and U = U + in satisfying (15.91) . (15.111) and solving the 
problems (15.36D - (15.401) . 


5.4. Problem (T2). Acting d* to the first equation in (15.391) and using the second equation, we 
derive that 

- A'ip p + H7ip p = 9o96(r p - r b ) + 9095 ip*dop + g&(U, U~, ip, DU~, DU, Dip, D 2 ip), (5.41) 

with go = /jo <75 + d*go and p 7 = — 9o96 < 0. Here A' is the standard Laplace operator on S 2 . Then 
using the third equation in (15..371) . we get 


- A\ip*p) + MV>*p) + M*P*dop) = 9s(U, U~,ip, DU, DU~,DiP, D 2 U, D 2 U~,D 2 iP , D 3 iP), (5-42) 

where pg = —909295 < 0 and g 8 = —A'g 2 + g 7 g 2 + 909295 + 92 d*go- We remind here again that 
DU and D 2 U represent first order and second order derivatives of ip*u,ip*p,ip*E,ip*A(s ) on S 2 . 

By the second equation in (|5.39l) . using the divergence theorem, and recall that f s2 ip p vo\ = 0, 
we have 

r p — r b = J^ (95 ip*(d 0 p) + 95 (U, U~,ip, DU, Dip)} vol. (5.43) 

Substituting this into the third equation in (15.371) . we then obtain 

ip p = — (ip*p-98 [ ip* (d 0 p) vol + g 7 (U, U~,ip,Dip)\ , (5.44) 

92 \ Js 2 J 

with // 8 = < 0 and g 7 = ^ J s2 g 5 vol - g 2 . 

We now formulate the following Problem (T2). 

Problem (T2): Find ip and U = U + U^ that solve (15.36[) . ()5.45[) . (|5.46|) . (|5.38j) . (j5.47[) 
and (15.401) . 


' L(p) = F(U,DU,D 2 p) in 

< p = pi -p£ on 

-A '( 4 >*p) + p 7 (ip*p) + 99 (ip*d 0 p ) 

= g 8 (U, U~,ip, DU, DU~ ,Dip, D 2 U, D 2 U~,D 2 iP, D 3 iP ); 

r p ~r b = -4^/52 (^95ip*(dop) + g 5 (U,U~, ip, DU, Dip)^ vol, 

< V = ^ {ip*P ~ 98 f s 2 ip*(doP) vol + g 7 (U, U~,ip, Dip)) , 
ip = ipP 4- r p ; 


(5.45) 


(5.46) 


d(ip*u’) = -±dg 0 (U,U , Dip), 

d*(ru') = 95 r(dop) + 90 V + 96 (r p ~ n) + g 5 (U, U~,ip, DU, Dip) 


Lemma 5.2. Problem (T2) is equivalent to Problem (Tl). 
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Proof. We only need to prove j s2 if p vol = 0, the third equation in (|5.37l) . and the first equation 
in (|5.39l) . Obviously the second equation in (15.471) and the definition of r p in (15.461) imply that 
f s 2 if p vol = 0. And the third equation in (|5.37l) . namely, (I5.27p . follows directly from (15.46[) . 

Now from the third equation in (|5.45l) . namely (|5.42p . and (15.271) . we get (15.411) . which means 
that, with the help of the second equation in (15.471) . 

d* (Vh/; p - no~ So) = 0. (5.48) 

The first equation in (15.471) says that the 1-form /loif* (vf) + go is closed. Since the first Betti number 
of S 2 vanishes, there is a function <f with f s2 </>vol = 0 so that d(f = poif*{u') + go- Furthermore, 
if f)' also satisfies dcf' = /^-oV 7 *(^- / ) + 9o, then d((f — (f r ) = 0, and we infer that <j>' = <f + k, with 
k € M. So from (I5.48p we have —A '{if p — (f) = d*d(if p — cf) = 0. By maximum principle, if p — cf 
is a constant. Since both of them have mean zero, we conclude that if p = cf. This shows that if p , 
hence ip, also solves the first equation in (15.391) . □ 

5.5. Problem (T3). The above equations and boundary conditions are formulated in At J. Sup¬ 
posing if € C 4,q ( 5 2 ), we introduce a C 4,Q -homeomorphism T : ( x°,x') € A4 4 ) i->- (y°,y ’) € H = 
(n,ri) x S 2 by 

y° = - -- 7 - 7 - 77(74 - n) + r b , y' = (y\y 2 ) = x’ = (x\x 2 ) (5.49) 

n — if{x') 

to normalize A4.J to 11. We set Ho = { p b} x S 2 and Hi = {rr} x S 2 . They are respectively the 
image of S^ and M\. Then 5H = Ho U Hi. We use i to denote the embedding of Ho in H. 

To avoid complication of notations, in the following we still write the unknowns in y-coordinates 
as (u,p, p, s, E) etc.; namely, we write (T~ 1 )*p as p, (T -1 )*^ as E, etc. and write 

y*u = ?1 r . b . fw°(^ -1 (y)) + -- —u a d a if(y')\ d y0 + u a {y))d ya (5.50) 

7 T - if(y) V n-r b ) 

still as u. (Note that (\H*u)“ = u a for a = 1,2. So there is no confusion.) By this convention if*E 
becomes i*E in y-coordinates. We have 

A ' x p = A' y p + 0(l)(D 2 pDif) + 0{l)(DpD 2 if) + 0(l){DpDif), 

Hence Problem (T2) could be rewritten as the following Problem (T3), where we use F or g to 
denote the higher order terms in the y-coordinates, and dj means 

Problem (T3): Find if G C 4, "(S' 2 ) and U = U+ U that solve the following problems 
(15.511) (15.561) . The initial data «q in (I5.56P is the vector field corresponding to the 1-form 
u' 0 on S 2 obtained from (15.551) . 


D U E = 0 in H, 

E = E~ — Eff on Hq; 


(5.51) 
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' L(p) = F(U,DU,D 2 p) in Q, 

p = pl-p b on ^ 1 , 

-A '(i*p) + p 7 (i*p) + p 9 (i*d 0 p) 

= g 8 (U, DU, Dip, D 2 U, D V, D 3 V>); 

r p -r b = - 4 ^ ,/ S 2 (55 i*(d 0 p) + g 5 (U, U~,ip, DU, Dif >)) vol, 

V = {i*P - P8 fs 2 i*(d 0 p ) vol + y 7 (U, U~,ip, Dtp)) , 

ip = ipP + r p ; 

f D u A(s) = 0 in fl, 

\i*(I(7)) = /i4 (^ p + r* - r 6 ) + 54(17, [/-, D^); 

fdK) = -^d5o(t/,C/-,^), 

\d*(u' 0 ) = P 5 i*(d 0 p) + p 6 p.> p + (r p — r h ) + g${U, U~, ip, DU, Dip) 

( G(D u u + , d a ) = 0 in Q, a = 1,2, 

lyin 0 = '“o- 

Remark 5.3. Recalling (I4.22p . in (|5.52|) we should have 

L (p) - ~j^y A 'p + (Ky 0 ) - l ) d lp + ^o b (t(y°)) d op + ^ 2 e ^y°))p 
+^^-d 1 (t(y°))E + ^^d 2 (i(y°))i(?) = F = F + F\, 


(5.52) 


on S 2 


(5.53) 

(5.54) 

(5.55) 

(5.56) 


( y °) 2 


(y 0 ) 2 


(5.57) 


where the coefficients are known functions of y°, and F = F(U, DU, D 2 p, Dip, D 2 ip) is the higher 
order term appeared below (14.221) . which is now writing in the y-coordinates; and 

A = 0{l){ip - r b )(D 2 p + Dp + U) + 0(l)(D 2 pDip) + 0{l)(DpD 2 ip) + 0{l){DpDiP). (5.58) 

We also note that 


98 = 98 + 0(ip - r b )Dp, g 7 = g 7 + 0(ip - r b )Dp, 

95 = 95 + 0(ip - r b )Dp, 54 = 94i 9o = go, 9e = Po95 + d*9o- 


(5.59) 

(5.60) 


Remark 5.4. Using (14.331) . we could obtain, for a = 1, the specific expression of (|5.56l) in a local 
coordinates in 12: 


W + 2ub y ) u 1 = -yyydip + fi(U, ip, Dip, Dp) in 12, 

on i2 0 . 


I u^d. 

\ u1 = u b 

Here u^dj is \k*u given by (|5.50|) . u' 0 = Uq d\ + u 2 d 2 , and 

/ 2 u b (y°) 2 u b (x° 


(5.61) 


fi(U, ip, Dip, Dp) = 


V y° 


1 , °1P 
u + 


P V(y 0 ) 2 (a : 0 ) 2 


1 


yO _ ri 


—diip doP + h^U). 


p(x 0 ) 2 yn — ip(y') 

By changing the local coordinates (y^yy 2 ), we know that u 2 solves a problem similar to (15.611) . 


(5.62) 
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5.6. Problem (T4). Since the elliptic problem (15.5211 is coupled with the hyperbolic problem 
(15.5411 . we need to further reformulate Problem (T3) equivalently as the following Problem (T4). 
We firstly note that by Lemma 15.21 Problem (T3) is equivalent to Problem (Tl), provided that 


if G C 4,a (S 2 ) and \if — r b \ < min 


n -r o ri — r b 


h^ 


4 4 

So we may replace the boundary condition in (|5.54|) by (15.301) . reads now 

i*(A(s)) = —i*(p ) + 34 - — 32 , 

M 2 M2 

where 32 = 32,34 = 34 , and the resultant problem is still equivalent to Problem (Tl). 
We now consider the Cauchy problem (|5.54jl : 


(5.63) 

(5.64) 


D u A(s ) = 0 

= + 34 -^32- 


in Q, 


(5.65) 


For the vector field u defined in Q, as in the proof of Theorem l6.ll we consider the non-autonomous 
vector field ^u(y°,y') defined for y' = (y-^y 2 ) G S 2 and y° G [r b ,ri\. For y G S 2 , we write the 

integral curve passing ( r b ,y ) as y / = (p(y°,y), which is a C k,a function in 11 if u G C k,a (Cl) and 

u° > 5 for a positive constant 5, and k G N. For fixed 3 0 , the map : S 2 —>• S' 2 , y y' = ip(y°,y) 
is a homeomorphism. Note that <p rb is the identity map on S 2 . 

Lemma 5.3. Suppose that u = u°d 0 + v! G C 0,1 and u° > 6. There is a positive constant 

C = C(5,r 1 — r b ) so that for any y' G S 2 and yo £ [n>, r i], it holds 

|(yyO'V - y'l < c |K|| c o ( q) • (5.66) 

Proof. Let y = 0/y>) _ V- Then |</y(y) -y\< f?° ^(s, <p(s, y)) ds < C ||w , |lc'0(o ) as desired. □ 
We write the unique solution to the linear transport equation (|5.65l) as follows: 

A{s){y) = A(s)(y°,y') = A(s){y°,<p y0 (y)) = A(s)(r b ,ip rb (y)) = A(s)(r b ,y) 

= (i*A(s))(y) = (i*4(s))((^o)-V)). (5.67) 

Hence, recall that the entropy is a constant behind the shock-front for the background solution, we 
have 

A(s)(y) = i*(A(s))(y') +^A(s)(r b ,((p y0 y 1 (y , ))-A(s)(r b ,y'^j 

= -i*{p) + 94 -— 92 +(A(s){r b ,(tp 0 )~ 1 (y')) - A(s)(r b ,y')Y (5.68) 

H2 Ll2 V / 


M2 


Set 


F 2 = - 


34 - ^-32 + ( A(s){r b , ( (fyO) 1 (y')) - A(s){r b , y')) | , 


(5.69) 


which is a higher order term (note that d a A(s ) itself is small, and <p y o is close to the identity map 
since u! is nearly zero, so \(<p y o )” 1 (y / ) — y'\ is small by (|5. 661) 1. Then we could write the elliptic 
equation (I5.57|) as 


1 A'p + (t(y°) - 1 )dlp + ^b(t(y°))d 0 p + 


( V °) 2 
Pb(y°) 
' (y °) 2 


y u 


(y 1 


e{t(y°))p + — Pb } y n 2 d 2 (t(y°))i*(p) 


0\2 


M 2 (y 1 


d 1 {t{y°))E + F + F 1 +F 2 . 


(5.70) 
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Now set 


ei(y°) 

= (V 0 ) 2 ( t(y°) - l) < 0, e 2 (y °) = 4 y°b(t(y 0 )), 

(5.71) 

e3(y°) 

= e(t(y 0 )), e 4 (y°) = —(p b (y°)yd 2 (t(y 0 )), 

P2 

(5.72) 

e 5 (y°) 

= -Pbiy^dytiy 0 )), 

(5.73) 

/ 

= f(U, DU , D 2 p, Dxl>, D 2 ^) = (y°) 2 (f + F 1 + F 2 ). 

(5.74) 

Then eauation (I5.70D simolv reads 


£(p) = -A 'p + ei(y°)dlp + e 2 (y°)d 0 p + e 3 (y°)p + e±(y°)i*(p) = e 5 (y°)E + /. 

(5.75) 


Note that there is a nonlocal term e 4 (y°)i*(p)- So problem ()5.52|) can be reformulated as follows: 


'm = e 5 (y 0 )E + f(U^,DU,D 2 p,D^,D 2 ^) in fi, 
p = p\ — p^ on 

-A \i*p) + (i*p) + n 9 (i*d 0 p) 

= gs(U,U-^ 1 DU,Di/ J ,D 2 U,D 2 ^,D 3 ^). 


We then state Problem (T4), which is equivalent to Problem (T3), as can be seen from the above 
derivations. 

Problem (T4): Find i/) € C 4 ’ a (S 2 ) and U = U^ + U defined in that solve problems 
(153T1) . (pTmD and (lQ3l) - (lT56l) . 


5.7. A linear second order nonlocal elliptic equation with Venttsel boundary condition. 

In this subsection we study the linear nonlocal elliptic equation (I5.75|) subjected to a Venttsel 
boundary condition on 

£(p) = / in fi, 

<p = hi(y') on fii, (5.77) 

^ -A \i*p) + PY{i*p) + p 9 {i*dop) = h 0 (y') on O 0 - 

Here / € C k ~ 2,a (£l),hi € C fc,Q (Hi) and h 9 € C k ~ 2 ’ a (Q o) are given nonhomogeneous terms, and 
k = 2, 3, • • • . This problem is formulated according to problem (15.761) . 

We remark that Venttsel condition of second order elliptic equations was proposed by A. D. 
Venttsel in 1959 in the studies of probability theory [25] • It is quite interesting to see that it 
appears naturally in the studies of transonic shocks. Later Y. Luo and N. S. Trudinger established 
the classical solvability for the linear and quasilinear Venttsel problems in a series of papers (see 
m for the linear case). However, our problem has two different characters comparing to these 
classical results. The first is that the elliptic operator £ in the domain H contains a nonclassical 
nonlocal term; the second is that the coefficients of the zero-th order terms may change signs. So 
we could not use directly the classical maximum principles or energy estimates, and we need the 
S-Condition to avoid some possible spectrums. 
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5.7.1. Uniqueness of solutions in Sobolev space i7 2 (H) n H 2 ( Ho). We firstly study under what 
conditions a strong solution p in Sobolev space i7 2 (H) with i*p £ H 2 (Ul q) to problem (15.771) is 
unique. To this end, we consider the homogeneous problem 

£(j5) = 0 in O, 

< p = 0 on Hi, (5.78) 

-A '(i*p) + /i 7 (i*p) + p 9 (i*d 0 p) = 0 on fi 0 . 

Remark 5.5 (Regularity). By Sobolev Embedding Theorem (Theorem 7.26 in [171 p.171]), H 2 (Q) 
is embedded into C°’^(Cl), so the solution is bounded; H 2 (Q 9 ) is embedded into C* 0,1 (Oo), so 
i*p £ By Trace Theorem, i*(d 9 p ) £ H 2 (fi 0 ). Now considering e 3 (y°)p + e A (y°)i*(p) 

as a nonhomogeneous term, the Schauder estimate (Corollary 6.7 in [T71 p. 100]) and a standard 
existence theorem of Dirichlet problems (like Theorem 6.13 in [IT]) imply that for any a £ (0, |), 
p belongs to C 2,a (Cl \ Ho). Then by Theorem 6.17 in [17], we infer that p £ C°°(Cl \ Ho) H C 0, z (Cl). 
We write the boundary equation in ([5.78j) as 

A'(i*p) = p 7 i*(p) + p 9 i*(d 0 p). (5.79) 

Note that the right-hand side belongs to iT^Ho). Applying the H s regularity theory (cf. Theorem 
11.1 in (21, Chapter 5, p.442]) and Sobolev Embedding Theorem, we have i*p £ Hz(£l o) CC 
C 1 ’? (Ho). Thus by intermediate Schauder estimate (cf. Notes of Chapter 6 in [T7] p. 138]), we have 
p £ C°°(£l U Hi) n C 1 ^ (fi). Hence the right-hand side of (|5.79|) lies in C 0, 2 (Ho)- By Theorem 6.14 
in m p. 107], as well as uniqueness of H 1 solutions of A' on S 2 (modulo a constant), we get that 
i*p £ C 2, 2 (Ho). Then Theorem 6.14 in [[171 p. 107] implies that p £ C 2 ' 2 (H), hence i*p £ C 3 ’2(Ho). 
Thanks to Theorem 6.19 in [TT] p.lll], we get that p £ C ,3 ’5(H), then i*p £ C 4 ’^(Ho), and hence 
p £ C 4 ’2(H), etc. Using this boot-strap argument, we deduce that if p £ H 2 (Ul) with i*p £ 77 2 (Ho) 
is a strong solution to ([5. 781) . then p is a classical solution and belongs to C'°°(H). 

Now we wish to show that p = 0. The idea is to use the method of separation of variables via 
the spherical harmonics. 

Let u n ^ m (y') be the eigenfunctions of —A' on S 2 with respect to the eigenvalues A n = n(n + 1) > 
0, n = 0,1,2, • • • , m = —n, ■ ■ ■ , — 1,0,1, • • • , n. (See Lemma A.4 in [7] p.2539] or j3[ Section 8.4, 
Theorem 1'].) Since {u n)m (y ')form a complete unit orthogonal basis of L 2 (S 2 ), we could write 

00 n 

P(V ) = ^2 V n,m(y°)Un,m(y'), (5-80) 

n= 0 m=—n 

with (recall that y = ( y°,y')) 

Vn,m(y°) = [ p(y°,y')u n , m (y')&S(y'). 

Js 2 

Here dS is the standard Lebesgue measure on S 2 . For p £ C k,a (Cl ) and k > 2, we easily deduce 
that u n ,m(y°) belongs to C k,a ([rb , h]) for all n, m , and it is also true that the series (I5.80P converges 
in C fc ^ 2 (H) (cf. [H Section 8.6, Theorem 3]). 

Substituting (15.801) into (15.781) . each v n ^ m solves the following nonlocal ordinary differential equa¬ 
tion: 


ei (y°)v",m + e 2 (y°)v' n:m + (e 3 (y°) + A n )v n ,m + e A (y°)v n , m (r b ) = 0, 
n = 0,1,2 ,--- ; m = -n, ••• ,n, y°£[r b ,r 1 \, 


(5.81) 
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and the two-point boundary conditions: 

v'n,m( r b ) + + Xn Vn t m(r b ) = 0, Vn.m^l) = 0. (5.82) 

hg 

We need to find conditions to guarantee that all v n ^ m are zero. 

Suppose that v' nm [r b ) = 0, then by uniqueness of solutions of Cauchy problem of ordinary 
differential equations, obviously one has that v n ^ m = 0 in [r&,ri]. 

If v nt m(r b ) / 0 for some n,m, we set w n = i c n ,m(y°) = v n , m (y°) / v n , m (r b ). Then it solves 

I ei (y°)< + e 2 (y°)w' n + ( e 3 (y°) + \ n )w n = -e 4 {y°), 

\w n (r b ) = 1, w' n (r b ) = w n (n) = 0. 

Definition 5.3. We say a background solution U b satisfies the S-Condition , if for each n = 
0,1, 2, • • • , problem (15.83[) does not have a classical solution. 


If the background solution U b satisfies the S-Condition, then all v n , m are zero, hence problem 
(I5.78P has only the trivial solution. Recall that a background solution is determined by the five 
parameters: 7 > l,r 6 G (r 0 ,ri),p s = p b (r b ) > 0 , p s = p b (r b ) > 0, M s = M fe + (r b ) € (0,1). We have 
shown in [7] that almost all background solutions satisfy the S-Condition. For example, we have 
the following lemma (c/. Lemma 2.6 in [7]). 

Lemma 5.4. For any given 7 > 1, M s € (0, 1), p s > 0 and p s > 0 , there is a set S C (ro,ri) 
of at most countable infinite points such that the background solution U b determined by 7, r b € 
(ro, r\) \ S, p s ,Ps, M s satisfies the S-Condition. 


Proof. Using ()3.4D and (|3.6I) . we note that all the coefficients ei,e 2 ,e 3 and e 4 could be solved, for 
the given constant 7 > 1 and initial data M s . p s . They are real analytic functions independent of 
r b . Now we change the independent variables y° to z given by z = - fj b . Then (I5.83H becomes 

V\ Vf) 


ei((n - r b )z + r b )w"(z) + (77 - r 6 )e 2 ((ri - r b )z + r b )w' n (z) 

+ (n - R,) 2 (e 3 ((ri - r b )z + r b ) + A n ^jw n (z) 

= -(n - r fe ) 2 e 4 ((ri - r b )z + r b ) z € [0,1], 

Wn (0) = l,w' n (0) = -(ri-r b )^, 
w n {l) = 0. 


(5.84) 


(5.85) 


Checking the expressions of ^7 and p 9 , they are analytic with respect to r b . So we see that the 
problem (|5.84l) depends continuously on r b S [r*o, r*i] and analytically on r b € [ro,ri). Hence the 
unique solution w n to this Cauchy problem is also real analytic with respect to the parameter r b . 
We write it as w n = w n (z ; r b ). Particulary, d n (r b ) = w n (l:r b ) is continuous for r b G [ro, n] and real 
analytic for r b € [ro, r 4 ). For given n = 0, 1,2, • • • , suppose now there are infinite numbers of r b so 
that 9 n (r b ) = 0. Then by compactness of [r Q , r 4 ], the function d n has a non-isolated zero point. So 
it must be identically zero and we have 6 n (r{) = 0. However, for r b = r\ , problem (|5.84[) is reduced 
to 

w"{z) =0, z € [0,1]; w n {0 ) = 1, ^(0) = 0. 

Hence w n ( 1) = 1, namely d n (ri) = 1 , contradicts to our conclusion that ’d n {fi) = 0. So for each 
fixed n, there are at most finite numbers of zeros of d n . Hence there are at most countable infinite 
numbers of r b so that the problem (|5.84l) and (|5.85l) may have a solution. The conclusion for (I5.83|) 
then follows. □ 
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5.7.2. Uniform a priori estimate in Holder spaces. By considering the nonlocal term e^y 0 )^(p) in 
£(p) as part of the non-homogenous term, and applying Theorem 1.5 in [TUI p.198] for the Venttsel 
problem (note that pg < 0) and Theorem 6.6 in [17] for the Dirichlet problem, with the aid of 
a standard higher regularity argument as in Theorem 6.19 of m. and interpolation inequalities 
(Lemma 6.35 in |T7t p.135]), we infer that any p G C k,a (Cl) (k = 2, 3) solves problem (15.771) should 
satisfy the estimate 

ll£llc fc ’ a (Q) — llpllc 0 (n) + ll/llc fc - 2 '“(ri) + ll^illc" c ’ Q (ni) + IIIIc ,fc - 2 ’“(f 2 0 ) )> (5.86) 

with C a constant depending only on the background solution Ub. 

Then, by (|5.86l) and a compactness argument, we have the a priori estimate: 

ll£llc fc ’“(Q) — c( ll/llc fc - 2 ’“(n) + IIMIc fc ’“(ni) + ||^o|lc" c - 2 '“(f 2 0 ) ) (5.87) 

for any C k,a solution of problem (I5.77|) . provided that the only solution to problem (15.781) is zero. 

In fact, to prove (|5.87|) , by setting K = ||/|| cfc - 2 ,a(n) + IIMc fe .“(L!i) + ll^-oIlc' fc - 2 .“(s^ 0 ) > we J ust 
need to show that there is a constant C so that 

IIpIIc 0 ^) < CK, (5.88) 

thanks to the inequality (I5.86p . Suppose that (15.881) is false. Then there are p^ G C k,a (Jl) so that 
||pW|| c °( 0 ) > riK n . Here K n is obtained from replacing / by /^ = £(p( n )), h\ by }i\ ' > = 

and kg by h\"' > = —A , (i*p( n ' ) ) + ^Hp^ + pgi* (dop^). By linearity of the problem, without loss 
of generality, we assume that for any n, llc°(n) = ^ Then K n -A 0 as n —>• oo and (15.861) 
implies that ||p^Ilc fc '“(n) — 2C. From Ascoli-Arzela Lemma, there is a subsequence {p^ nj ^}j of 
so that converges to a p^ in C k (yi) (k = 2,3). Thus by uniqueness of solution of 
problem (|5.77l) . taking j —> oo, we must have p$ = 0, contradicting to the fact that ||^llo°(f2) = 
limj 

—^OO Thus we proved (I5.88j) . 

5.7.3. Uniform a priori estimate in Sobolev spaces. Suppose now that p G 7f 2 (fl) D H 2 (Q o), which 
means that p G i7 2 (H) and the trace of p on Ho, namely i*p , belongs to H 2 {Ul o). Obviously our 
assumptions on problem (|5.77l) guarantee that / € L 2 (H), h\ € H 2 (S 2 ), and ho G L 2 {S 2 ). Then 
by Trace Theorem and Interpolation Inequalities of Sobolev functions (c/. Theorem 7.28 in i 17. 
p.173]), we have 

IK*T|Il 2 (S 2 ) ^ C I|p|l^i(f2) < e l|p|lH 2 (f2) + C( e ) IIpIIz,2(0) , (5.89) 

\\i*(d 0 p)\\ L2{s2) < \\i*(d 0 m H ^ s2) <C\\d 0 p\\ H s m <C\\p\\ H r m 

< £ \\p\\H 2 (n) + C( £ ) ll^lli 2 (n) > V e G (0,1). (5.90) 

Now applying Theorem 8.12 in m p.186] to the boundary equation in (15.771) . we have 

ll**Pll.ff 2 (S 2 ) — II**pIIl 2( S 2 ) + ||h 0 || L 2( S 2 ) + ||**(<9 op)|Il 2 (5 2 ) ) 

< IIpII_h 2 (q) + C(e) ||p|| L 2(Q) + C 11/1011^2(52) • (5.91) 

Using the same theorem to problem (15.771) . with given Dirichlet data i*p, it follows that 

\\p\\H 2 (n) + I Idlin' 2 (s 2 ) - c( II**£IIh 2 (5 2 ) + ll^i 11 h 2 (s 2 ) + IIpIIl 2 (q) + ll/II.L 2 (n)) 

< C £ \\p\\H 2 (fl) + C{£) ||p|| L 2 ( n ) + c(^ ||ho|| L 2 ( 5 2 ) + ||hl|| H 2 ( S 2 ) + ||/||x, 2 (Q) 
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Taking e = 1/(2 C), we get 

\\p\\H 2 (n) + W^PWmis 2 ) — IIpIIl 2 (G) + II^o||l 2 ( 52 ) + HMI# 2 ^ 2 ) + \\f\\L 2 (n) ) • (5.92) 

With an argument similar to the proof of (15.871) (c/. Lemma 9.17 in [171 p.242]), we deduce the a 
priori estimate 

\\p\\H 2 (n) + W*p\\h 2 {s 2 ) — \\^o\\l 2 (s 2 ) + ll^illiL 2 (s 2 ) + ll/llz, 2 (n) ) > (5.93) 

provided that the S-Condition holds. Here the constant C depends only on the background solution 
U b and r h ,r\. 


5.7.4. Approximate solutions. We now use spherical harmonic expansion to establish a family of 
approximate solutions to problem (|5.77D . 

For simplicity, we take h\ = 0 in the sequel. There is no loss of generality, since this accounts 
we replace p by p — h\, and / by / — 2h\, ho by ho + A'h\ — yjhi in problem (|5.77l) . 

We also set {f^}k to be a sequence of C 00 ^) functions that converges to / in C k 2, “(H), and 
C C°°(S 2 ) converges to ho in C k ~ 2,a (S 2 ). Now for fixed k, we consider problem (15.771) . 
with / there replaced by f^ k \ and ho replaced by h$ . 

Suppose that 

oo n oo n 

f (k \y) = Y^ fh%(y°) U n,m{y'), h ( 0 k) (y') = Y^ (4 fc) )n,mMn,m(y')- ( 5 - 94 ) 

n=0m=—n n=0 m=—n 

Then for p given by (|5.80p . each v njrn (y°) should solve the following two-point boundary value 
problem of an ordinary differential equation containing a nonlocal term: 


L n ,m(Vn,m) = <, m (y°) + P(y°H,m(y°) + Qn(y°)Vn,m(y°) 

= r{y 0 )v nt m(r b ) + f n ,m(y°) y° € [r b ,n], 

'Un,m(j'b) T 0“nVn,miXb) — ^ n,mi 

V n ^m{ri) = 0. 

Here we define 

e 2 (y°) , 0 \ e 3 (y°) + X n 0 e 4 (y°) ~ 0 fn%{y°) 


P(y ) = —To)’ qn(yV ) = 

d(y u ) 


: ( o\ _ ) f /O', _ 

ei (y u ) e 1 (y u ) e 1 (y u ) 


O'77. - 


fJ>7 + A 


n i _ 

i n nm — 


(h( k h 

Wo )n,r. 


P9 M9 

We will show that this problem is uniquely solvable. 

Let m be the unique solution of the Cauchy problem 

^/n,m(^n,m) — 0; — 1; m^Tb) — 9) 


and </? 2 m the unique solution to 


h'n,m{'^n,m) 0, 'Cn,m(cfc) 9 ■ '^n,mi.'^b^ 1* 


(5.95) 


Then by standard theory of linear ordinary differential equations, a general solution to problem 


I J n,m(.'Un,m) — c(y T /n,m(?/ ) 2/ € 

V n ,m(r b ) = Ci, v n,m( r b) = C 2 


(5.96) 
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is given by 


Vn,m(y ) = ciip n)7n (y ) + c 2 y n)m {y )+ / K(y ,s)(cir(s) + / n , m (s))ds, 


>r b 


(5.97) 


where 


K(t,s)± 




We claim that there are constants c\ and c 2 so that 
C 2 Q"nC 1 — hn,mi 

i) + Sri K(n,s)r(s ) ds) Cl + <#U(ri)c 2 = - JT K(n, s)f n , m (s) ds. 


(5.98) 


Actually, this is a linear algebraic system and we know that, under the S-Condition, the homo¬ 
geneous system has only the trivial solution. So there is one and only one pair (ci,c 2 ) solves 
(15.98)1 . Hence ((5.97)1 gives the unique solution to problem (15.951) . Note that f n . m G C°°([rb, rr]) as 
f(k) g C°°(Cl), and the coefficients p,q n ,r are all real analytic, so the solution v n ^ m (y°) belongs to 

Now for N € N, we define 

N n 


Pn(v ) — /* ' 'y ^ v n,m(y ) u n,m{y )j 


n =0 m=—n 


N n 


N n 


fN ) (y) = '52 X fn%(y°) u n,m(y'), = X X 

77 .— 0 777——71 77—0 777——77 

Apparently pat, € C 00 ^), and (/iq^at € C'°°(5 2 ). It is also easy to check that pat solves the 
following problem: 

£(Pn) = f n ] in 

Pat = 0 on f2i, (5.99) 

-A'(i*p N ) + y 7 {i*p N ) + p 9 (i*d 0 p N ) = (h^) N on £l 0 . 


5.7.5. Existence. By the estimate (15.93)1 . for any N, M € N with N < M, there holds 


WPm ~ PNWn^^ + W^iPM ~ Pn)\\ H 2^ o) < C 




L 2 (S 2 ) 


+ 


f(^) f(k) 

J M J N 


L 2 (t2) 


Recall that {h^)N —t h^ in L 2 (S 2 ) and —> /( fc ) in L 2 (H), we infer that {pat} (respectively 

i*Pn) is a Cauchy sequence in H 2 (Q) (respectively H 2 (£l q)). So there is a p^ G H 2 (£l) (respectively 
g(fc) g h 2 ( Qo)) an d Pat —>■ p^ in H 2 (Q) (respectively i*pN —>■ g( fc ) in i7 2 (17 q)) as IV —>• oo. By 
continuity of trace operator, we conclude that = i*p^ k \ Taking the limit IV —» oo in problem 
(15.990 . one sees that p( fc i is a H 2 (Q) n H 2 (Qq) solution to problem (15.770 . with / there replaced 
by f^ k \ ho replaced by h^\ Then by the same arguments as in Remark 15.51 p( fc i G C' 00 (I7) and of 
course it satisfies the estimate (15.870 . 

Now for the approximate solutions {p^}k, we use the estimate ((5.87)1 to infer that 

|p( fe ) 

< \\f\\c k - 2 ’ a ((i) + ll^o|lc fc - 2 ’“(f2o) )• 
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Hence by Ascoli-Arzela Lemma, there is a subsequence of { p that converges to some p € C k,a (£i) 
in the norm of C fc (fi). Taking limit with respect to this subsequence in the boundary value problems 
of p k , we easily see that p is a classical solution to problem (15.7711 . Therefore, we proved the following 
lemma. 

Lemma 5.5. Suppose that the S-Condition holds. Then problem (|5.77l) has one and only one 
solution in C k,ol (fil), and it satisfies the estimate (15.8711 . 

5.8. Solvability of Problem (T4). We now use Banach fixed point theorem to solve the transonic 
shock problem (T4), provided that the background solution U b satisfies the S-Condition. 


5.8.1. The iteration sets. Let ctq be a positive constant to be specified later, and 

= jt/> G C 4,a (S 2 ) : \\ip - r b || C 4, Q(S 2) < cr < £j 0 | 

be the set of possible shock-front. For any given if € /C CT , its position r p and profile if p also satisfy 

\r p - r b | < cr, ||^ p || C 4,«(s2) < 2cr - 

We write the set of possible variations of the subsonic flows as 

Os^{u = (E,p,s,u') : ||^|| 3 + ||**^|| 0 3 ,a (sa) <6<6 0 }, 

with 5o a constant to be chosen. The norm ||-|| fc appeared here is defined by (I4.35p . with A4 there 
replaced by fi. 

Given U~ satisfying ()5.15|) . for any if € K, a and U € Os, we construct a mapping 

T : K a X O s ->■ K a X Os, W,U)^$,U) 

as follows. One should note that a fixed point of this mapping is a solution to Problem (T4). 


5.8.2. Construction of iteration mapping. For any if € JC a and U € O 5 , we set 


U = U + lV 


-{y ~ n) + fi{y'),y’). 


r\ - r b 

Then with the known supersonic flow U ~, we could specify all the higher-order terms / and g bj 
appeared in Problem (T4). 

Bernoulli constant E. We firstly solve the linear problem (c/. (15. 511) ) 


(5.100) 


J D U E = 0 in 0, 
yE = E~ — E b on Ho- 

Note here that E~ is determined by U~, and E b is a constant. So by (|5.15j) . it holds 

||i*(£' — E b ) 11^3,= || {E — E b )\ s ^ || C 3, a ^ 2 ) < Cq£. (5.101) 

Hence we could easily get the unique existence of E € C 2,a {fii) (note that u 6 C 2 ,Q (H)) with 


E 


!, a /m - W l *( E - C ° £ - 


(5.102) 


c 2 .“(a) 

The constant Cq appeared here and below depends only on the background solution and r b ,r\. 
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Pressure p. We then consider the following problem on p (cf.. (15.761) 1: 

ft(p) = e 5 (y°)E + f(U,^DU,D 2 p,D^,D 2 ft in fi, 

< P = Pi-pt on ^i) 

-A'(i*p) + p 7 (i*p) + p 9 (i*d 0 p) 

= gs(U,U~^,DU,D^,D 2 U,Dft,Dft). 


(5.103) 


Here the non-homogeneous terms / and gs are determined by if E IC a and U = U + U b , with 
U E Os, and E is solved from (15.1001) . Then, since we assumed that the S-Condition holds, by 
Lemma [531 we could solve uniquely one p E C 3,a (Q) and it satisfies the following estimate: 


IIpII C 3 ’ a (Q) - Co[ WfWc^m + II58|IcL«(5 2 ) + 




*(s 2 ) 


E 


c 1 ' a (n) 


+ lh Pb llc 3 '“(ni) 


(5.104) 


Checking the definitions of / and gs, the right-hand side is finite; actually we have (see Lemma [6.11 
in the Appendix) 


ll/llc 1 ’ Q (f2) < 0(6" + + e 2 + e), 1 1 <7811 o 1 (S' 2 ) — E{5 2 + e 2 + + e). (5.105) 

So combining (15.1021) . (15.81) and (15.1041) . (15.1051) . one infers that 


llp|lc 3 ' Q (n) — Co (<5 2 + + £ 2 + • 


(5.106) 


Update shock-front if. With the specified higher-order terms g b and 57 , and p solved from 
([5.1031) . we now set (cf. (15.531) 1 

r p -r b = j (ji 5 i*(dop) + <?5 (U, U~,if, DU, Dif)j vol, 

‘ ft = — (i*p- p s [ i*(d 0 p) vol+ g 7 (U, U~,if,Dif)\ , 

P 2 V Js 2 J 

if = ft + f p . 

It follows easily that (using ([5.1061) 1 

ft 


(5.107) 


C°(S 2 ) 


+ \r p — r b \ < C 0 l 1135licks' 2 ) + ll 37 |lco(sa) + 


< C 0 (^ + cT 2 + e 2 + e). 


(5.108) 


For the C^ ,a estimate of ft , we note that i*p solves the third equation in ()5. 103[) . hence ft solves 
the following elliptic equation on S 2 (cf. ([5.41)1 ): 

- A'ft + p 7 ft = /r 0 M6 (r p - r b ) + p 0 p 5 i*8 0 p + g 6 (U, U~,if, DU~,DU, Dif, D 2 if). (5.109) 

Standard Schauder estimates on Chapter 6 ] yield that 


ft 


04.0(52) - C 0 (^ 11^11(70(52) + ll56|lc 2 .“(S 2 ) + IIpIIc 3 -“(5 2 ) + \ pP ~ r b 

< C'o(<5 2 + cr 2 + e 2 + e). 


One then has 


4>-r b 


C 4 -“(S 2 ) 


<C 0 [5 2 + a 2 + £ 2 + £ 


(5.110) 

(5.111) 
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We also need to show that 


[ $ p vol = 0. 

Js 2 


(5.112) 


In fact, integrating (15.1091) on S' 2 , and recall g$ = Kofi's + d*goi using Divergence Theorem and 
definition of r p — r b in (15.1071) . we have directly (I5.112p . 

Entropy H(s). Now we solve (c/. (15.541) . and note that H(s)j~) is constant) 

j D u A(s) = 0 in D, 

1 i*(A(s)) = /i 4 ( i>-r b ) +g i (U,U~,'iJ),Dil>) 


(5.113) 


to obtain the unique solution 7l(s). It also holds 


i*A(s 

A(s 


C 3 ’ a (S 2 ) 


<C 0 


- r b 


C 3 -“(S 2 ) 


+ II 54 II 


C 3 ’ a {S 2 ) 


_ <Co 
C 2 ’“(G) 


i*A{s) 


C 2 .“(S 2 ) 


< ColS^ + ^ + ^ + e). 


< Co (<5 2 + cr _ + £ 2 + e j ,(5.114) 


(5.115) 


Tangential velocity field u'q on Do. Next we study the tangential velocity u' {) on Do by (c/. 
(1535D ) 


on S 2 . (5.116) 


\d(u' 0 ) = -—dg 0 (U,U 

\d*(«o) = Ms i*(doP ) + Me V + Me (r p - r 6 ) + g 5 (U, U~,ip, DU, Dif}) 

By (15.1121) and the first equation in (15.1071) . we see the requirements in Theorem 16.21 are fulfilled. 
So one can solve a unique u' 0 on S 2 , and the following estimate is valid: 


I u' 


o|lc 3 '“(s 2 ) — c 'o(^||Mllc 3 .«(a) + 


i>-r b 


C 2 .“(S 2 ) 


+ ll5 , o|lc 3 .“(S 2 ) + Hs , 5|lc 2 .“(S 2 ) 


(5.117) 


< C 0 (<5 2 + u 2 + e 2 + e). 

Note that u' 0 is a 1-form on S 2 . Let u' 0 = u^dp be the associated vector field on S 2 . Then according 
to (15.1171) . Uq = Uq satishes the following inequality 

2 


£ 

0=1 


-.0 


C3.“(S 2 ) 


< C 0 (5 2 + u 2 + e 2 + e). 


(5.118) 


Tangential velocity u' in D. Finally, we solve the tangential velocity v! in D through (c/. (1 5.61 1) 1 

,ja s.i 1 o^bvy > ?. 1 _ 1 


WjU 1 + 2'^jj^-u 1 = —fijrpdip + fi{U,'ip,D'i/j) in D, 


lA = Un 


on Dq. 


(5.119) 


Here the Cauchy data Uq on Do is solved from (15.1161) . 

We obtain a unique u l in D by Theorem l6.il From (15.1041) and (15.1181) . it holds that 


7 ) 1 _ < 

1 “ llc 2 .“(n) — 


Co( IImIIc 3 -“(0) + ||/i| 


C 2 .“(f2) 


+ 


c 2 ’“(n 0 ) 


< C 0 (<5 2 + cr 2 + e 2 + e) . (5.120) 

By changing the coordinates system, we get u 2 and it also satisfies an estimate like ()5.1201) . 
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if — r b 

+ 

u 

+ 

i*U 

<c( 


C 4 ’“(5 2 ) 


3 


C 3 -“(S 2 ) V 


Conclusion. From the above six steps, we get uniquely one pair (U, if) and it follows from (|5. 101 D . 

flEUSD , dEM}, (EUD, <nrm . iEnm ietttbii and ' m that 

(5.121) 

Here C is a constant depending only on the background solution and r&, r\ . Now we choose C* = 4 C 
and Eq < min jl/(16(7 2 ), 1, ^/( 8 C')|. Then, for <5 = a = C*£, we have C^ 5 2 + a 2 + e 2 + e'j < 

5, \/e € (0, £o)) and the estimate (|5.121l) shows that U G Oc,e and if G K,c„e- Hence we construct 
the desired mapping T on /Cc* e x Oc.e- 

5.8.3. Contraction of iteration mapping. What left is to show that the mapping 

T : Kc.e X Oc,e ->• /C ae x 0 C .e, 0, &) ^ (Vb U) 

is a contraction in the sense that 


flV - ip) 

1 


< 


+ 

C 3 ’“(5 2 ) 

-y>( 2 ) 


[/(!) _ t>(2) 

+ 

C 3 >“(S 2 ) 


+ 


i*0W - t> (2) ) 


17 (1) - C/ (2) 


+ 


C 2 .“(S 2 ) 

T(C/ (1) - C/ (2) ) 


C 2 .“(S 2 ) 


a 1 


= t 2 Q; ( 5 - 122 ) 


provided that £q is further small (depending only on the background solution). Here for j = 1,2, 
and any if^ G /Cc* e , G Oc, £ , we have defined U^) = T(if^\ U^). 

To prove (|5.122|) . we set if = if^ — and U = t/l 1 ) — tj( 2 \ For k = 1,2, we also use the 
notations 

(U~)M = U ~, T(f/+)W = t/+|^ w , 

(C/+)( fc ) = ( U+){ ri ~^ ){ - y '\ y° - r 6 ) + if W (y'),y'), m + ([/+)«. 

ri rjj 

By (|5.15l) and analyticity of U±, the mean value theorem implies that there is a constant C 
depending only on the background solution so that 


z*(£/+) (1) -i*(t/+) (2) 

W)«-W)«|| cl „ (il) 

(U- - U b -)M - (U- - ( 7 t -) (2) 


C fc -“(S 2 ) 

< C 


< C 


- ipw 

ifW _ ^( 2 ) 

< Ce 


C k ’ a (S 2 ) 


C 2 -“(5 2 ) 

i. We note that E solves the following problem (c/. (15.1001) ) 

(D uW E + D uW _ u(2) E^ = 0 in n, 
{E = E-\ s ^-E-\s^ on n 0 
By mean value theorem, (|2.1D and (|5.15l) . there holds 


C k ’ a (S 2 ) 
ifW - if® 


k = 1,2,3,4, 
C 2 -“(S 2 ) ’ 


(5.123) 

(5.124) 

(5.125) 


i*E 


C 2 ’ a (S 2 ) 


<9n E 


0 ^ II C 2 ’ a (M) 


ifW - v> (2) 


< C 
Ce 


(u?) dgE- 

if^ _ ^(2) 


C 2 ’ a {M) 

C 2 ’ a {S 2 ) 


C 2 ’ a (S 2 ) 

ifW _ ^(2) 


C 2 ’ a (S 2 ) 


(5.126) 
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Then using Theorem 16.11 and (|5.124l) , one has 


E 

< 

i*E 

+ 



£-(2) 


CL“(G) 


CL“(S 2 ) 


CL“(G) 



< Ce 

< Ce 


</;(!) -^(2) 

^ 


C 2 ’ a (S 2 ) 


+ CC,e 


c 2 ’ a {n) 
[/(!) _ #(2) 


+ 


uw _ f/(2) 


cL“(n) 

< CeQ. 


+ 


U+ {1) - U+ 


( 2 ) 


C 2 ’ a (S 2 ) 

Step 2. Next we seek an estimate of p, which solves (c/. (|5. 1U3I) 1 

£(p) = e 5 (y°)E + / (1) - / (2) in fi, 


CL“(fi) 

(5.127) 


p = 0 


on fli, 


(5.128) 


-A'(i*p) + p 7 (i*p) + p 9 (i*d 0 p) = 5g 1} - 9s ) - 


Here for k = 1,2, 

f {k) = f(U {k \^ {k \DU {k \D 2 p { - k \D^ k \D 2 'il: { - k) ), 

g {k) = g 8 {U W ,{U-) {k \^ k \DU W ,D^ k) ,D 2 U W ,D 2 ^ k \D 3 ^ k) ) 

By Lemma [5751 Lemma [6721 and (I5.127p . direct computation yields 

E _ + / (1) - / (2) _ + 


\\p\\c^m ^ c v 

Step 3. From (15.1071) . we see that 


C“(H) 


-(i) -(2) 

98 -38 


C a (S 2 ) 


< CeQ. (5.129) 


4 ye 7 sl v« i 'W ) + ^ 1, -^ )) ) TOl ’ 

V = j- 2 (i*P - P8 f s 2 i*(d 0 p ) vol + (<^ 1} - 4 2) )) 
^ + fP. 


Here, for k = 1,2, 


# =g 5 (U {k \(U-)^,^ k \DU^,D^), 

=Mu {k \{u-i k \^ k \D^). 

Then we have the following estimate via (15.1291) , Lemma 16.21 and some straightforward computa¬ 
tions: 


V c(s 2 ) + n c l ^cnn) + 


-( 1 ) -( 2 ) 

% 


C(VL) 


+ 


-( 1 ) -( 2 ) 
97 - g 7 


C(n) 


< CeQ. 

By (15.1091) . note that ip p also solves 

- A'^ p + pr4) p = p 0 p 6 f p + g 0 g 5 i*d 0 p + g^ - gf ^, 

with 


(5.130) 


_(fc) 

% ~ " 


= g 6 (U^ k) ,(U-^ k \^ k \(DU-Y k) ,DU^ k \D^ k \D 2 ^), k = 1,2, 
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it follows that, from (15.1301) and (15.1291) . 




C 3 ’ a (S 2 ) 


< C 


i/fi 


C(S 2 ) 


+ \r p \ + l|p|lc 2 ’“(a) + 


-( 1 ) -( 2 ) 
96 -96 


C 1 ’“(S 2 ) 


< CeQ. 

This and (15.1301) imply that 

Step 4■ From (I5.113p . one has 




C 3 .“(S 2 ) 


< CeQ. 


( 2 ) 


(5.131) 


D uW A ( s ) + D uW _ u (2)A(s) =0 in Q 

i*(A(s )) = 4> + gP ~ 94\ 


where 


g {k) = g 4 (U {k \ {U~p\p k \ DpV). k = 1,2. 
By (15.1311) and Lemma 16.21 we have 


It is also easy to show that 


C 2 '“(S 2 ) 


< C 




C 2 .“(s 2 ) 


+ 


-( 1 ) -( 2 ) 
94 -9\ 


C 2 ’ a (S 2 ) 


< CeQ. 


Ap) 

< 

i*(A(s)) 

+ 

uW-u® 


1 ) 

s 




C^ts 2 ) 


C^iP) 



C 2 ’ a (fl) 


(5.132) 

< CeQ. (5.133) 


Step 5. Next we turn to (15.1161) to find that the difference of tangential velocity field on f2o solves 
Jd(«o) = -9o ] ), 


(d*(h' 0 ) = g 5 i*(d 0 p) + gei/j + (g 5 ~ 9 5 ) 

where, for k = 1,2, 

§P = Hu {k \{u-p\p k \D^ k) l 
§P = g 5 (U^ k \(U-p\^ k \DU^,D^). 

By Theorem 16.21 we have the estimate 


on S 2 , 


-( 1 ) -( 2 ) 
00 -00 


+ I|p|lc 2 ^(t2) + 




+ 


-( 1 ) -( 2 ) 
05 -05 


Po||c?2,«(5 2 ) - C y -u -u C 2,« (S 2) ■ oo ■ ' CT'TS 2 ) ' " u (^’“(S 2 ), 

< CeQ. (5.134) 

Step 6. From (15.1191) . u l , the first component of the difference of the tangential velocity in 12, solves 
(uQ^djii 1 + 2 u ^o - u 1 + ((uQA') — (uQ^^Jdjtu^) 1 

= ~ pwlyoydiP + j^yi p( p w p p P drP^ + (fP - fP) in 12, 


U 1 = Un 


on No¬ 


where 


fP = h(U^,^ k \D^), k = 1,2. 
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So there holds 

'HI 


u ^c 1 > c *(Q) — C 


Ic 2 ’“(f2) 


+ 


/l (1) - /l (2) 


_ + C*£ 
CL“( a) 


J7 (1) - C/ (2) 


CL“(G) 


< CeQ. 

There is a similar estimate for the second component, namely u 2 . 


(5.135) 


Conclusion. Now summing up the inequalities (15. 1261)(15.1271) (15. 1291) (15. 131 II (l5.132D(j5.133D(l5.134|) 
and (j5. 1351) . we get 


$ 

+ 

U 

+ 

i*U 


C3.“(S 2 ) 


2 



C 2 >“(S 2 ) 


< C'eQ, 


which implies (15.1221) if e € (0,£o) and C'eq < 1/2. Finally, by Banach fixed point theorem, we 
infer Problem (T4), hence Problem (T), has one and only one solution in K-c*e x Oc t£ ■ This finishes 
the proof of Theorem 15.11 


6. Appendix 

We provide here some results used in this paper, together with some details on the estimates of 
higher-order terms. 

6.1. Solvability and estimate of transport equations. We consider the following Cauchy 
problem of a transport equation for the unknown E in M = (ro,7"i) x M , where M is a smooth 
closed surface: 

D U E + aE = f in M., E = Eq on Mq = {ro} x M. (6.1) 

The main result is: 

Theorem 6.1. For a fixed number a € (0,1) and k G N, suppose that the vector field u = u °<do + u' 
and the functions a, f belong to the Holder space C k,a (M), and Eq € C k,a (M), and furthermore, 
u° has a positive lower bound 5 in A4. Then there is uniquely one solution E to problem (EH), and 
there is a positive constant C = C(M, \r\ — ro|,<5, ||u|| C fc, Q , ||a|| C fc, a ) so that 

\\-^\\c k ’ a (M) — ll-^oll C k ’ a (M) + \\f\\c k ’ a (M) ) • (6-2) 

Proof. 1. Since u° > S > 0 in A4, we may rewrite problem (16.11) as 

\i! CL f 

d 0 E + —rrE + -jrE = -77 in M, E = E 0 on M 0 . 

■u u n u 

Hence thee is no loss of generality by assuming that u° = 1 in the sequel. 

2. For u' = u a d a , consider the Cauchy problem of ordinary differential equations: 

<\ T ' 

— = u'{t,x'), x\ro) = x € M. 
at 

So x'{t) is the integral curve of the (time-dependent) vector field u' on M, which passes through 
the point x on M when t = ro- Since k > 1, and v! is bounded, by theorem of ordinary differential 
equations, there is one and only one such solution for t € [ro,ri], and the solution depends on x 
with the same regularity as v! depending on x. (See, for example, [26l Section 13].) We may also 
write the solution as x' = (p(t, x) = (ft(x). The transform pt '■ M M, x i-a x' and its inverse both 
belong to C k ’ a (M). 
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3. For the transport equation 8qE + u'E + a{x)E = f(x ) of unknown E(x) (we write here that 
x = ( t,x '); that is, x° = t ), since E{x) = E{t,x ') = E(t,ipt(x)), we have ^ + a(t, (x))E = 
f(t,ip t (x)). The solution is given by 

E(x) = E(t,ip t (x)) = exp ( - / a(s,^ s (x))dsj 

' Jro ' 

x( [ E 0 {x) + exp a(r,</? r (x))dT^/(s,</? s (x))ds^ . (6.3) 

This proves the existence. 

4. Since products and compositions of C k,a functions are still in C k,a , and note that x = 

(ip t )~ l {x') € C k,a (M), (</? t ) _1 is C k+1,a with respect to t. we deduce that the right-hand side of 
(16.31) belongs to C k,a ( A4), and the estimate (16.21) follows. □ 


6.2. Solvability and estimate of div-curl system on sphere S 2 . The following result is a 
refined version of that appeared in [7], so we omit the proof here. 

Theorem 6.2. For a € (0,1) and k € N, there exists a unique 1-form u in C k+l {S 2 ) that solves 


= x, d*cu = if, 

if X is a 2-form in C k,a (S 2 ) and if is a function in C k,a (S 2 ), and there hold 


(6.4) 



(6.5) 


Furthermore, there is constant C > 0 so that 

ll w llc fc + 1 -“(5 2 ) — c( llxllc fc ’“(S 2 ) + \\' l l’\\c k ’ a (S 2 ) )• 


( 6 . 6 ) 


6.3. Estimates of higher-order terms. We list below the estimates of some higher-order terms, 
which were used in Section 15.8.21 and Section 15.8.31 in studying the iteration mapping T■ 


Lemma 6.1. There is a constant C depending only on the background solution Ub so that if 
if € lC a ,U € Os, and U~ satisfying (15.151) . then one has the following inequalities: 


IMIc3,a(S2) 

ll/llcL“(n) 

ll-^ll C 2 -“(f2) 

llffo|lc 3 . a (S 2 ) 


< c[5 2 + a 2 + e 2 + e 

< C(5 2 + cr“ + £ 2 + e), 

< C^d 2 + a 2 + £ 2 + £ 
C (5 2 + a 2 + £ 2 + £ 


||ff5|lc 2 .“(S 2 ) — C (V + <J~ + E 2 + £^j . 

< C(6 2 +(T 2 +£ 2 +£^, 


Wa&Wctx 


{so ^ 


||57|lc°(S 2 ) — C(j) 2 + (T~ + £ 2 + £^J . 

< c(d 2 + £ 2 + a 2 + £ y 


11^8 llc'i.c 


(so ^ 


(6.7) 

( 6 . 8 ) 
(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 


In the expressions of these higher order terms, we have set U = U + (y° — rf) + if(y'),y'). 










44 


LI LIU, GANG XU, AND HAIRONG YUAN 


Proof. 1. We firstly consider the terms g k (i.e., g k ) defined through (|5.26I) - (I5.29I) for k = 1,2,3,4. 
Checking (2.27) in [3 p.2522] shows that g k are of the form 

g k = 0(l)i*{U)(iP-r b ) + 0(l)i*{U) 2 + 0(l)i*{U)r(U~ -U>) 

+0(l)ijj*(U~ - U b ) + 0(l)(Dip)i* (U) + 0(1)(V> - r b f. (6.15) 

It follows easily that 

||<7fc|lc' 3 ’ Q! (S 2 ) — + S 2 + eS + £ + 5ct + o' 2 ) < C{^5 2 + <r" + e 2 + £^j. (6.16) 

Thus we proved (16.71) . 

2. The term / is defined by (15.7411 . For F 2 (see (|5. 691) 1. we note that 

A( s )( r b, (p y o)-\y')) - 4l(s)(r fe , y') = 0(l)(D(i*U))U, (6.17) 

hence according to the proof of Lemma 15.31 its C' 1, “(f)) norm is controlled by C5 2 . Then using 
(16.151) . we find that 


F 2 


C L“(f2) 

For the term F\ (see (I5.58p l. it can be written as 


<C{5 2 Fa 2 + e 2 + e). 


(6.18) 


F\ = 0(l)(iP-r b )(D 2 p + Dp + U) + 0(l)(D 2 ^)(Dp) 

+0(l)(D^)(D 2 p) + 0(l)(Dfj)(Dp). (6.19) 


Hence 



< Cda < C(5 2 + a 2 ). 

cL“(n) 


( 6 . 20 ) 


Next we consider the term F (see (14.221) 1. For F 2 (see (14.271) 1. taking into account of the transform 
T (see ([5.491) 1. it is of the form 


F 2 = 0(1)(U ) 2 + 0(1)(DU)U + 0(l)(D 2 p)U + 0(1)(DU) 2 , (6.21) 

where the quantities 0(1) may contain — r b , Dpi or D 2 ip. Therefore one has 

ll^llcv^) < C<5 2 . (6.22) 

For H 3 (see (14.141) 1. it can be written as H 3 = 0\U 2 , with Oi a bounded quantity containing 
D 2 p, (if - r b ),Di/> and Dp. Then F\ = Oi(H 2 + U(DU) + (DU) 2 ) and 

Pill &,«(&)< C 6 2 . (6.23) 

Summing up the above estimates ()6.181)(16.201)(16.221)(16.23p . we get (16.81) . 

3. The term f\ is given by (15.621) . where h\ is defined by (14.431) . We may write 


fi = 0(l)U 2 + 0{l)U(rJ> - r b ) + 0(l)(Dp)(^ - r b ) + 0(l)(Dp)(Di/>), 


where 0(1) may contain Dip and ip — r b . Therefore it is straightforward to check that 
which implies (16.91) . 

4. The terms g k for k = 0,5, 6 , 7 ,8 were defined in ([5.591) and (15.601) . We firstly consider go (see 
(15.24p 1. It is of the form 

go = 0 (l)(i*U ) 2 + 0(l)(FtJ)r(U~ - U b ~) + 0 (l)ip*(U~ - U b ), 
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where 0(1) may contain Dp). Therefore 

ll3ollc3,«(s 2 ) — + £ )- 

Thus we proved (16.1011 . 

5. Next we consider 55 , where 55 was given by (15.351) . and G was defined by (15.331) . We see 

55 = 0(l)(i*U)(D(i*U)) + 0(l)i*U(P)-r b ) + 0(l)Dfi(i*U) 2 
+0{l)Dp){i*U) + 0(l)Dfi(i*U)(D(i*U)) + 0(l)(i*f7 ) 3 
+0(l){i*Uf{D(i*U)) + 0{l){i*U) 2 + 0(l)Di/)(fi*U) 2 (D{i*U)) 

+0{l)(fi) - r b ) 2 + 0(l)pj*(U~ - f7 fe “) + 0(1)52 + 0 ( 1)54 + 0(l)(p> - r b )Dp. 

Hence we have 

||55|lo2.«(5 2 ) < C(S 2 + cr" + e 2 + e) 

as desired. 

6 . Since ge = /ro 5 s + d* 5 o, the estimate (| 6 . 12 l) follows directly from (16.101)(16.111) . 

7. We note that g 7 = f s2 g 5 vol - 52 , so 

97 = 0(1)55 + 0 ( 1)52 + 0(1)(V> - r b )(Di*p) (6.24) 

and from (16.111) . (16. 16|) . we get (16.131) . 

8 . Finally we consider 5 §, where gs = —A'g 2 + 5752 + 9296 - So (16.141) follows from (| 6 . 12 l) and 
( 1061 ) . □ 

Lemma 6 . 2 . Under the assumptions stated in Section [5.8.3l there is a constant C depending only 
on the background solution so that 


sf-S? 

C 2 ’ a (S 2 ) 

< 

CeQ , 

(6.25) 

/ (1) - / (2 

C“(0) 

< 

CeQ, 

(6.26) 

/i (1) - /i (2) 

CL“( n) 

< 

CeQ, 

(6.27) 

-(1) -(2) 
So -So 

C 2 ’ a (S 2 ) 

< 

CeQ, 

(6.28) 

-(1) -(2) 
95 -95 

C^iS 2 ) 

< 

CeQ, 

(6.29) 

-(1) -(2) 
So - So 


< 

CeQ, 

(6.30) 

_(1) _(2 
9) ~ 9) 

) 

C(S 2 ) 

< 

CeQ, 

(6.31) 

-(1) -(2) 
98 -98 

C“(S 2 ) 

< 

CeQ. 

(6.32) 


Recall that Q is defined by (15.1221) . 

Proof. To prove (16.251) . we use (16.15f) . Returning to (2.27) in [7], some tedious computations show 
that for each coefficient 0(1), it holds ||0(1)( 1 ' ) — 0(1)( 2 ) || C 2 , a (g 2 ) < CQ. Hence, as an example, 
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since 


0(l) (1) (^ (1) )*(C/ _ - U~) - 0(l) (2) (^ (2) )*(C/~ - U b ) 

= 0(1)W (V 1] W - U-) - (v> (2) )*(u - - u>j) 
+(o(i)« - o(i)^)(^nu- - u b ~), 


by (15.1251) we have 


0(l) (1) (^ (1) )*(C/“ - U b ) - 0(1) (2) (V> (2) )*([/" 




C 2 ’“(S 2 ) 


< CeQ. 


The other example is to consider 

0(l) (1) r(17 (1) )(^ (1) - r b ) - 0(l) (2) r(17 (2) )(^ (2) - r b ) 

= 0(l) (1) r([/ (1) )(^ (1) - V> (2) ) + (V> (2) - r b )0(l) w i*(U (1) - C/ (2) ) 

+(t/> (2) - r 6 )(0(l) (1) - 0(l) (2) )i*[/ (2) , 


and we have straightforwardly that 

0(l) (1) i*(l7 (1) )(V> (1) - r b ) - 0(l^ 2) i*(tJ^)(^ 2) - r b ) 


C 2 ’ a (S 2 ) 


< CeQ. 


Since all other terms are of similar form, they could be estimated similarly and hence we omit the 
details. □ 


Acknowledgments. Li Liu was supported by National Natural Science Foundation of China 
(NNSFC) under Grants No. 11101264, No. 11371250. Gang Xu was supported by NNSFC under 
Grants No.11101190, No.11371189 and No.11271164. Hairong Yuan (the corresponding author) 
was supported by NNSFC under Grants No.11371141. Yuan would also like to thank the Wuhan 
Institute of Physics and Mathematics, Chinese Academy of Science for the hospitality, where the 
manuscript was finished. 


References 

[1] Bae, Myoungjean; Feldman, Mikhail. Transonic shocks in multidimensional divergent nozzles. Arch. Ration. 
Mech. Anal. 201 (2011), no. 3, 777-840. 

[2] Benzoni-Gavage, S.; Serre, D. Multidimensional Hyperbolic Partial Differential Equations: First-order Systems 
and Applications. Oxford Mathematical Monographs. Clarendon Press, Oxford, 2007. 

[3] Bleecker, D.; Csordas, G. Basic Partial Differential Equations. International Press, Boston, 1996. 

[4] Chen, Chao; Xie, Chunjing. Three dimensional steady subsonic Euler flows in bounded nozzles. J. Differential 
Equations 256 (2014), no. 11, 3684-3708. 

[5] Chen, Gui-Qiang; Chen, Jun; Feldman, Mikhail. Transonic shocks and free boundary problems for the full Euler 
equations in infinite nozzles. J. Math. Pures Appl. (9) 88 (2007), no. 2, 191-218. 

[6] Chen, Gui-Qiang G.; Deng, Xuemei; Xiang, Wei. Global steady subsonic flows through infinitely long nozzles 
for the full Euler equations. SIAM J. Math. Anal. 44 (2012), no. 4, 2888-2919. 

[7] Chen, Gui-Qiang G.; Yuan, Hairong. Local uniqueness of steady spherical transonic shock-fronts for the three- 
dimensional full Euler equations. Commun. Pure Appl. Anal. 12 (2013), no. 6, 2515-2542. 

[8] Chen, Shuxing. Stability of transonic shock fronts in two-dimensional Euler systems. Trans. Amer. Math. Soc. 
357 (2005), no. 1, 287-308 (electronic). 

[9] Chen, Shuxing. Stability of a Mach configuration. Comm. Pure Appl. Math. 59 (2006), no. 1, 1 -35. 

[10] Chen, Shuxing. Transonic shocks in 3-D compressible flow passing a duct with a general section for Euler systems. 
Trans. Amer. Math. Soc. 360 (2008), no. 10, 5265-5289. 

[11] Chen, Shuxing; Yuan, Hairong. Transonic shocks in compressible flow passing a duct for three-dimensional Euler 
systems. Arch. Ration. Mech. Anal. 187 (2008), no. 3, 523-556. 












THREE-DIMENSIONAL STEADY COMPRESSIBLE EULER SYSTEM 


47 


[12] Courant, R.; Friedrichs, K. O. Supersonic flow and shock waves. Applied Mathematical Sciences, Vol. 21. 
Springer-Verlag, New York-Heidelberg, 1976. 

[13] Dafermos, C. M. Hyperbolic Conservation Laws in Continuum Physics. Third edition. Grundlehren der mathe- 
matischen Wissenschaften, vol. 325. Springer-Verlag, Berlin Heidelberg, 2010. 

[14] Dzhuraev, A. Systems of equations of composite type. Pitman Monographs and Surveys in Pure and Applied 
Mathematics, 44. Harlow: Longman Scientific & Technical; New York: John Wiley & Sons, Inc., 1989. 

[15] Fang, Beixiang; Liu, Li; Yuan, Hairong. Global uniqueness of transonic shocks in two-dimensional steady com¬ 
pressible Euler flows. Arch. Ration. Mech. Anal. 207 (2013), no. 1, 317-345. 

[16] Frankel, T. The Geometry of Physicist, An Introduction, 2nd Ed. Cambridge University Press, Cambridge, 2004. 

[17] Gilbarg, David; Trudinger, Neil S. Elliptic partial differential equations of second order. Reprint of the 1998 
edition. Classics in Mathematics. Springer-Verlag, Berlin, 2001. 

[18] Li, Jun; Xin, Zhouping; Yin, Huicheng. Transonic shocks for the full compressible Euler system in a general 
two-dimensional de Laval nozzle. Arch. Ration. Mech. Anal. 207 (2013), no. 2, 533-581. 

[19] Luo, Yousong; Trudinger, Neil S. Linear second order elliptic equations with Venttsel boundary conditions. Proc. 
Roy. Soc. Edinburgh Sect. A 118 (1991), no. 3-4, 193-207. 

[20] Liu, Li. Unique subsonic compressible potential flows in three-dimensional ducts. Discrete Contin. Dyn. Syst. 27 
(2010), no. 1, 357-368. 

[21] Liu, Li; Yuan, Hairong. Stability of cylindrical transonic shocks for the two-dimensional steady compressible 
Euler system. .1. Hyperbolic Differ. Equ. 5 (2008), no. 2, 347-379. 

[22] Liu, Li; Yuan, Hairong. Uniqueness of symmetric steady subsonic flows in infinitely long divergent nozzles. Z. 
Angew. Math. Phys. 62 (2011), no. 4, 641-647. 

[23] Liu, Li; Yuan, Hairong. Steady subsonic potential flows through infinite multi-dimensional largely-open nozzles. 
Calc. Var. Partial Differential Equations 49 (2014), no. 1-2, 1-36. 

[24] Taylor, M. Partial Differential Equations I: Basic Theory. 2nd edition. Springer-Verlag, New York, 2011. 

[25] Venttsel, A. D. On boundary conditions for multi-dimensional diffusion processes. Theor. Probability Appl. 4 
(1959), 164-177. 

[26] Walter, Wolfgang. Ordinary differential equations. Graduate Texts in Mathematics, 182. Readings in Mathe¬ 
matics. Springer-Verlag, New York, 1998. 

[27] Wang, Ya-Guang; Yuan, Hairong. Weak Stability of Transonic Contact Discontinuities in Three-Dimensional 
Steady Non-isentropic Compressible Euler Flows Z. Angew. Math. Phys. (2014), on line. (DOI: 10.1007/s00033- 
014-0404-y) 

[28] Weng, Shangkun. A new formulation for the 3-D Euler equations with an application to subsonic flows in a 
cylinder. (2012) Arxiv: 1212.1635v2. 

[29] Yuan, Hairong. On transonic shocks in two-dimensional variable-area ducts for steady Euler system. SIAM J. 
Math. Anal. 38 (2006), no. 4, 1343-1370 (electronic). 

[30] Yuan, Hairong. Transonic shocks for steady Euler flows with cylindrical symmetry. Nonlinear Anal. 66 (2007), 
no. 8, 1853-1878. 

[31] Yuan, Hairong. A remark on determination of transonic shocks in divergent nozzles for steady compressible Euler 
flows. Nonlinear Anal. Real World Appl. 9 (2008), no. 2, 316-325. 

[32] Yuan, Hairong. Examples of steady subsonic flows in a convergent-divergent approximate nozzle. J. Differential 
Equations 244 (2008), no. 7, 1675-1691. 

L. Liu, Shanghai University of International Business and Economics, Shanghai 201620, China 

E-mail address: llbaihe@gmail.com 

G. Xu, Faculty of Science, Jiangsu University, Zhenjiang, Jiangsu 212013, China 

E-mail address: gxu@ujs.edu.cn 

H. Yuan (corresponding author), Department of Mathematics, East China Normal University, 

Shanghai 200241, China 

E-mail address : hryuan@math. ecnu .edu.cn 


